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ABSTRACT 

These lecture notes are addressed to PhD student and/or researchers who want a general 
overview of microscopic approaches based on mean-field and applied to nuclear dynamics. Our 
goal is to provide a good description of low energy heavy-ion collisions. We present both formal 
aspects and practical applications of the time-dependent Hartree-Fock (TDHF) theory. The TDHF 
approach gives a mean field dynamics of the system under the assumption that particles evolve 
independently in their self-consistent average field. As an example, we study the fusion of both 
spherical and deformed nuclei with TDHF. We also focus on nucleon transfer which may occur 
between nuclei below the barrier. These studies allow us to specify the range of applications of 
TDHF in one hand, and, on the other hand, its intrinsic limitations: absence of tunneling below 
the Coulomb barrier, missing dissipative effects and/or quantum fluctuations. Time-dependent 
mean-field theories should be improved to properly account for these effects. Several approaches, 
generically named "beyond TDHF" are presented which account for instance for pairing and/or 
direct nucleon- nucleon collisions. Finally we discuss recent progresses in exact ab-initio methods 
based on the stochastic mean- field concept. 



* Lecture given at the "Joliot Curie" school, Maubuisson, September 17-22, 2007. A french version is available at 
the URL http://www.cenbg.in2p3.fr/heberge/EcoleJoliotCurie/coursannee/cours/CoursSimenel.pdf . 
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I. INTRODUCTION 
A. General considerations 

Heavy-Ion accelerators have given important information on nuclear reactions with stable nu- 
clei. More and more precise measurements shed light on the interplay between reaction mechanisms 
and the internal structure of the two reaction partners. This interplay is perfectly illustrated by 
low energy reactions like fusion. Fo r insta, nce, fusion cross sections are influenced by vibrational 
and rotational modes of the nuclei Das98l | . New radioactive low energy beams facilities such as 



SPIRAL2, open new opportunities for these studies. Theoretically, microscopic models which incor- 
porate both dynamical effects and nuclear structure in a common formalism should be developed. 
Dynamical theories based on the mean-field concept are of the best candidates. 

The present lecture notes give a summary of actual progress in mean-field transport models 
dedicated to Heavy-Ion reactions at low energy. The starting point of different approaches, i.e. 
Time-Dependent Hartree-Fock with effective interaction, is described extensively in section |TI1 
while examples of applications are given in section lllli TDHF should be extended to describe the 
richness of phenomena occurring in nuclear systems. In section ITV] we introduce transport theories 
that incorporate effects such as pairing correlations and/or direct nucleon-nucleon collisions. 

This lecture can be read at different levels. The reader interested in applications (section IIIip 
can skip most of the formal aspects described in section |TI1 Section IIVI is relatively technical and 
requires a good understanding of the formal aspects described in section [Hi Finally, minimal notion 
of second quantization which are used in these notes are summarized in appendix 1X1 

B. Microscopic Non-relativistic approaches 

Typical examples of quantum microscopic theories used to describe static properties of atomic 
nuclei are: various Shell Models, mean-field and beyond theo ries (H artree-Fock-Bogoliubov, 
Random-Phase- Approximation, Generator-Coordinate- Method...) [RinSOl ] or algebraic approaches 



like the Interactive-Boson- Model Iac91 |. 



Why should we also develop quantum microscopic approaches to describe nuclear reactions? 

• By itself, the N-body dynamical problem is a challenging subject. 

• Most of the information one could get from nuclei are deduced from nuclear reactions. There- 
fore a good understanding of the reaction is mandatory. 

• As already mentioned, it is necessary to develop a common formalism for both static and 
dynamical properties. 

• Microscopic theories have only few adjustable parameters (essentially the eff^ective interac- 
tion), and we do expect that the predicting power is accordingly increased compared to more 
phenomenological approaches. 

• On opposite to most of macroscopic approaches which are specifically dedicated to a given 
mechanism, all dynamical effects should be "a priori" included. For instance, TDHF can be 
used either to study reactions between two nuclei or collective motion in a single nucleus. 

• Using microscopic theories, we also expect to be able to deduce dynamical information on 
the behaviour of nucleons in the nuclear medium like for instance in-medium nucleon-nucleon 
cross sections 

Last, microscopic approaches that will be described here are non relativistic. Relativistic effects 
are not expected to affect atomic nuclei because the typical velocities of nucleons are much less 
than the light speed (w/c)^ ~ 1/10. 
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C. What means independent particles? 

It is relatively common to use words like "mean- field theories" or "independent particles sys- 
tems" as well as their complement "beyond mean- field" or "correlated states". Let us start with 
a discussion on this terminology. Then, we will be able to define the starting point of most of 
microscopic approaches dedicated to the description of nuclei. 



1. Independent and correlated particles 

We first assume two identical free fermions, one in a state \^) and the other one in a state [ly). 
Except for the Pauli principle which forbids = {u), the fact that one fermion occupies |//) is 
independent from the fact that the other occupies |z^). We say that the two particles are "indepen- 
dent". The state vector associated to the two particles system accounts for anti-symmetrization 
and reads = (|1 : ^, 2 : i/) — |1 : i^, 2 : /\/2. Such a state is indifferently called independent 
particle state or Slater determinant (see appendixE]). The concept of independence could be easily 
generalized to particles. A system will be called here "independent particle system" if it could 
be written in a specific basis as an anti-symmetric product of single-particle states. In second 
quantization form (appendix [Aj), such a state writes 




On the opposite, we call "correlated" a state that cannot be written in terms of a single Slater 
determinant, i.e. \jp) = Ca \(j)a) where each Cq, gives the probability of each configuration. 
Let us for instance consider a state describing two particles that decomposes onto two Slater 
determinants, i.e. {tp >= ci|/iiz/i > +C2\fJ-2i^2 >■ Then the concept of correlation becomes obvious 
because the occupation of one specific state affects the occupation probability of the other state. 
In this particular example: 

• If one particle is in the state then the other is in 

• If one particle is in the state |i^2)) then the other is in |/i2)- 



2. Mean-field approximation 

The aim of nuclear mean-field theories is to describe self-bound nuclei in their intrinsic frame 
where wave- functions are localized (on opposite to the laboratory frame). A possible description 
of a self-bound localized system in terms of Slater determinants could be constructed from single- 
particle wave-functions of an harmonic oscillator or a Woods-Saxon potential. These potentials are 
then interpreted as effective average mean-fields that simulate the interaction between particles. In 
other words " Each nucleon freely evolves in a mean- field generated by the surrounding nucleons" . 
It is worth mentioning the necessity to consider intrinsic frame is very specific to self-bound systems 
like nuclei. For instance, for electronic systems, electrons are always considered in the laboratory 
frame since they are automatically localized due to the presence of atoms and/or external fields. 

In this lecture, we consider more elaborated mean-field, generally "called" self-consistent mean- 
field, like those found in Hartree-Fock (HE) and Hartree-Eock-Bogoliubov (HEB) theories. In the 
first case, the mean-field is calculated from occupied single-particles while the second case is more 
elaborated and requires the notion of quasi-particles (see section [IV Ap . In section [II B| TDHE 
equations are obtained by neglecting correlations. Last, although it is not the subject of the 
present lecture, it is worth mentioning that the introduction of effective interactions for nuclear 
mea n-field theories leads to a discussion on correlation much more complex than presented here 



(see |Sto07l | for a recent review). 
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3. Mean-free path and justification of mean-field approaches in nuclear physics 

Can independent particles approximation give a good description of nuclear systems including 
reactions between two nuclei? The justification of mean-field theories is largely based on the 
empirical observation that many properties vary smoothly for nuclei along the nuclear charts: 
single-particle densities, energies... Based on this consideration, macroscopic and mean-field models 
have been introduced and turns out to be very successful to describe nuclei. Therefore, by itself, 
the predicting power of independent particle approximation justify a posteriori its introduction. 

The large mean-free path of a nucleon in nuclear matter (larger th an the size of the nucleus 
itself) gives another justification of the independent particle hypothesis Boh69l | . This implies that 



a nucleon rarely encounter direct nucleon-nucleon collision and can be, in a good approximation, 
considered as free. This might appear surprising in view of the strong interaction between nucleons 
but could be understood as a medium effect due to Pauli principle. Indeed, the phase-space 
accessible to nucleon after a direct nucleon-nucleon collision inside the nucleus is considerably 
reduced due to the presence of other surrounding nucleons (essentially all states below the Fermi 
momentum are occupied). However, if the relative kinetic energies of two nucleons increases, 
which happens when the beam energy increases, the Pauli principle become less effective to block 
such a collision and the independent particle approximation breaks down. 

4- Symmetries and correlations 

Small remarks against the intuition: independent particle states used for nuclear systems contain 
correlations. This could be traced back to the fact that some symmetries of the original Hamiltonian 
are generally broken. For instance, nucleons described within the nuclear mean-field approach are 
spatially correlated because they are localized in space. This is possible because mean-field is 
introduced in the intrinsic frame and translational invariance is explicitely broken. Indeed, if we 
do not break this symmetry, then the only mean-field solution would be a constant potential and 
associated wave- functions would be plane waves. We then come back to the free particle problem 
which are not self-bound anymore. 

We illustrate here an important technique which consists in breaking explicitly symmetries to 
incorpor ate corr elations which could hardly be grasped in an independent particle picture (see for 
instance jBen03l |). Among the most standard symmetries explicitly broken generally to describe 
nuclear structure, we can quote breaking of rotational invariance which authorizes to have deformed 
nuclei and help to recover some long range correlations^. Gauge invariance (associated to particle 
number conservation) is also explicitly broken in HFB theories in order to include short range 
correlations like pairing (see section IIV Ap . The latter approach will still be called " mean- field" . 
It however goes beyond independent particle approximations by considering more general Many- 
body states formed of product of independent quasi-particles (a summary of terminology and 
approximations considered here is given in table HI). 

Last, it is worth mentioning that any brok en syinmetry should normally be restored. This is 
generally done by using projection techniques [RinSOl . lBla86l . [Ben03l | . 



5. Theories beyond mean-field 

In nuclear physics, mean-field is often considered as the "zero" order microscopic approximation. 
Many extensions are possible (several "beyond mean- field" approximations will be presented in 
section lIVp . These extensions are generally useful to include correlations that are neglected at the 



^ Let us consider an elongated nucleus ("cigare" shape), the fact that one nucleon is at one side of the cigare implies 
necessarily that other nucleons should be on the other side. This should be seen as long range correlations that 
affects nucleons as a whole. 
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mean- field level. As we will see, "beyond mean-field" approximation are absolutely necessary to 
describe the richness of effects in nuclear structure as well as in nuclear dynamics. In particular, not 
all correlations could be incorporated by only breaking symmetries and often one has to consider 
the state of the system as a superposition of many independent (quasi) particles states. 

The previous discussion clearly points out that mean-field approaches might include correlations. 
Nevertheless, the terminology that is generally used (and that we continue to use here) is that 
non- correlated state will be reserved to Slater determinant states. A correlated state then refers 
to a superposition of Slater determinants. Table [I] summarizes different approaches that will be 
discussed in this lecture and associated type of correlations included in each approach. 



Name 


Approximation 


Variational space 


Associated observables 


TDHF 


mean- field 


indep. part 


one-body 


TDHF-Bogoliubov 


m.-f. + pairing 


indep. quasipart. 


generalized one-body 


Extended-TDHF 


m.-f. + collision (dissipation) 


correlated states 


one-body 


Stochastic-TDHF 


m.-f. + collision 
(dissipation-f fluctuations) 


correlated states 


one-body 


Time Dependent Density Matrix 


cm. + two-body correlations 


correlated states 


one- and two-body 


Stochastic Mean Field 
(Functional integrals) 


Exact (within statistical errors) 


correlated states 


all 



Table I: Summary of microscopic approaches presented in these notes 



D. Effective interaction and Energy Density Functional (EDF) 

We will use in the following a rather standard approach to the nuclear many-body problem. 
Starting from a microscopic two-body Hamiltonian and using second quantization, we introduce 
the Hartree-Fock theory. It is however worth mentioning that the use of most recent realistic 
two-body (and normally three-body) interactions will not lead to reasonable results (if any) at 
the Hartree-Fock level. In view of this difficulty, the introduction of effective interactions adjusted 
to nuclear properties was a major break-down. These interactions are not directly connected to 
the original bare interaction and are expected to include effects (in particular in-medium effects) 
that are much more involved than in the standard Hartree-Fock theory. In that sense, mean-field 
theories in nuclear physics have many common aspects with Density Functional Theories (DFT) in 
condensed matter. However, in nuclear physics we often keep the concept of effective interactions. 

The most widely used interactions are contact interactions (Skyrme type) and finite range 
interactions (Gogny type). The second type of interactions is still too demanding numerically to 
perform time-dependent calculations and only Skyrme forces are nowadays used for TDHF. To 
relax some of the constraints due to the use of effective interactions, the more general concept of 
Energy Density Functional (EDF) is introduced. In that case, the static and dynamical properties 
of the s ystem are direc tly obtained by minimizing a functional of the one-body density, denoted 
by E[p] [Dob07l . lBmOi l. 
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E. The N-body problem: basic formalism 



The evolution of a Many-Body state \tp) is given by the time-dependent Schrodinger equation 

ih^\i;)=H\i;). (2) 

This equation can equivalently be formulated in term of a variational principle. In that case, 
the quantum action is minimized between two times to and ti with respect to the wave-function 
variations 



dt (VI H - in- 
to ^'^ 



0. (3) 



The Many-body Hamiltonian decomposes into a kinetic term and a two-body interaction term 
(for the sake of simplicity, we will not consider higher order interactions) 

" = Hi.f)- (4) 

i=l i>j=l 

In second quantization (appendix |A]), the Hamiltonian writes 

H = ^ tij alcLj + ^ Vijkl a\a\aiak (5) 

ij ijkl 

where matrix elements associated to the kinetic energy and anti-symmetric two-body interaction 
are given respectively by 

tij = ^ i^\p'\j) (6) 

Vijkl = Vijkl — Vijlk (7) 
Vijkl = {l:i,2:j\ v{l,2) \l:k,2:l). (8) 

Note that here, the two-body state \1 : i,2 : j) is not anti-symmetric. This notation means that 
the particle "1" is in the state \i), while the particle "2" is in the state \j) (for further details on 
notations see appendix |A|) . 

Let us again summarize the overall goal: Even if we were able to solve the exact Many-Body 
problem (Eq. ([2])), which is not the case anyway in nuclear physics, we do not need to have all 
the information contained in the exact wave-function to understand physical processes related 
mainly to one-body degrees of freedom, like collective motion or fusion reactions. Construction of 
microscopic models are therefore guided by two principles: 

• we focus exclusively on what we really need, 

• we do relevant approximations for the considered problem. 

A strategy of approximation consists in focusing on specific degrees of freedom. This is generally 
equivalent to minimize the variation al prin ciple ([3|) on a restricted subspace of the total Hilbert 
space of Many-body wave functions [BlaSSl ] . The selection of this subspace is crucial and is driven 
by the physical process considered. 

Among the microscopic transport theories, the TDHF meth od, proposed b y Dirac in 1930 Dir3Cll | 



as an extension of the Hartree-Fock static mean- field theory [Har28l . lFoc30l ] is a tool of choice. It 



corresponds to a mean-field theory where the only input is the effective two-body interaction. In 
practice, TDHF Equations can be derived by restricting the variational space to Slater determi- 
nants. 

In the following, formal and practical aspects of TDHF are presented. Then examples of ap- 
plications to reactions close to the fusion barrier are given. Successes and limitations of TDHF 
will illustrate physical effects missing in this approximation. "Beyond mean-field" theories that 
incorporate these effects are finally discussed. 



9 



II. DYNAMICAL MEAN-FIELD THEORIES 



Let us recall the two basic questions: 

• What is the relevant information for the description of nuclear reactions ? 

• What approximations can (and cannot) be done to treat this information ? 

The knowledge of nuclei trajectories during collisions, their shapes or particle numbers already 
gives a good understanding of the reaction mechanisms. All these quantities are related to one- 
body degrees of freedom (see appendix [B]). Therefore, if we are able to give a realistic description 
of the one-body density matrix, we will also correctly reproduce these observables. As we will 
see, the reduction of the information to one-body degrees of freedom is intimately connected to 
the independent particle approximation discussed in the introduction. This will answer the second 
question. Starting from the Schrodinger equation, the equation of motion for the one-body degrees 
of freedom (contained in the evolution of the one-body density) is obtained. Then connections 
with the independent particle approximation are discussed. Finally, we describe practical aspects 
related to the resolution of TDHF. 

Readers not interested in formal details can jump to Eq. ()34p which gives the most general form 
of the one-body density evolution for any Many-Body system. 



A. Expectation values of one-body observables 

1. General expression 

Let us consider a system of N particles IV') (eventually correlated) and a one-body operator 
F = X]j=i /(O (see Eq. ()Bip ). Starting from general expression of a Many-Body wave-function 
(Eq. ()A22p ) and using the closure relation (Eq. (|A28P ). we obtain 



1 r ^ 

= ^ jdii..AiMdi[..Ai'^ riii-^N) m-^N) E^^i-^^i /(^) 1^1-^^) (9) 



i=l 



where spin and isospin quantum numbers are included in the notation ^ = (rsr). Using the Wick's 
theorem (see appendix IA3p . leads to 



Nl 



m'l) 



(eii/ieo (eiieUi) 



Nl 



iiNWi) ■■■ {iN\iu mm {iNW,+i) ■■■ {^nw^) 



(10) 



Using the anti-symmetry of each state of the basis, it is possible to reduce, in integral Q, the 
action of / to the first label "1" only. Noting /(^,^') = {i\f\^') the matrix elements of /, we 
obtain 



(11) 
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2. One-body density matrix 

We introduce the one-body density matrix p^^^ associated to the state IV')- Its matrix elements 
in the basis {^} are given by 

p^'H^,o = N jd^i-d^N r{e^2..-CN) ipm-^N). (12) 

Accordingly, equation ([TT]) can be simply expressed as 

(F)^ = yded^' P^'He,0 = Tr[p(i)/]. (13) 

Therefore, all the information required to estimate a one-body observable is contained in the one- 
body density matrix. Such a matrix can always be associated to any system of N particles^. In the 
following, we will essentially consider properties of the one-body density and omit exponent 
Using second quantization, matrix elements of p are defined as (see Eq. ()C2p ) 

Pij = {tl^l a^-cii 14)) = {a^jCLi)^. (14) 

Eq. dUD can indeed be recovered using (jATSl) . (|A22]) and ([A28l) 

p(c,e') = {a\e)mh 

= ^ JdCi...dCNdC[...de^ ri^i-^N) HC'i-C'n) {^i-Cn\ a^(e')«(0 WvC'n) 

= yd6...dC7vd^^...de^v n^'a-^N) V'(ce2-e^v) (6-e7vie2-e^> 

= N ld^2:.d^N r{^'^2-4N) H^a-^N). (15) 

In the following, the one-body density components are introduced in a specific single-particle 
basis. For instance (see section Hi Pp . TDHF calculations are often performed in coordinate repre- 
sentation. In that case, using equations ()A6p and (|A7p . we have 



p(r.r,r'/r') = (V^| at(r'/r') a(r.r) |^) = Y.{ala,)^ ^t'V) (^f (r) (16) 



where single-particle wave- functions tpf^ are defined in Eq. ()Aip . 

Note that, we can always write the one-b ody den sity matrix as an operator acting in the Hilbert 
space of the single-particle wave-functions^ Rin80l | 



P = T.P^i (17) 



In particular, this operator depends on the system wave- function {ip) and then it is a funtion of 
the time during a collision. 



It is worth mentioning that the discussion presented here can be generahzed to M-body observables. In that 
case, we can define the M-body density that contains all the information on Af-body (and below) quantities (see 
appendix [Cj. 

^ This definition of operators p should not be confused with the definition of one-body operators given in appendix [B1 
since it is not defined in the space of many-body wave function. 
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Let us come back to one-body observables. The second quantization simplifies the calculation 
of their expectation values. Using Eq. (|B3p . we have directly 

(F)^ = (V^l ^ a\ a, IV^) = ^ p,, = Tr (p/) (18) 

which is nothing but Eq. ()13p . Therefore we see that the evolution of any one-body observable 
can be obtained from the evolution of the one-body density matrix. 

3. The independent particle case 

Let us consider the case of a system described by a Slater determinant \4>) = '^^^ 
associated one-body density is denoted by p. In appendix [El we show that any Slater can be 
considered as a vacuum (called HF vacuum) for specific operators written as linear combination 
of the and a. We can thus apply the Wick's theorem. Methodology associated to the Wick's 
theorem as well as contractions are given in appendix [Fl This technique (Eq. (jFSp ) implies that 
only occupied states contribute to the summation in Eq. ()16p . We finally end with 

N 

p{rsT, v's'r') = ^i:'V) <(r). (19) 
n=l 

Therefore, for Slater determinants, the knowledge of occupied states is equivalent to the knowledge 
of one-body density matrix. The specificity of independent particles systems is that all the infor- 
mation is contained in the one-body density. This is nicely illustrated by the fact that any M-body 
density matrix could be expressed as an anti-symmetric product of one-body densities leading to 
vanishing M-body correlation matrices at all orders (see appendix [C]). 

For dynamical problems where the Many-Body state is assumed to stay in an independent 
particle state at all times (this is the case of TDHF presented in section III Bp , it is sufficient to 
only follow the one-body density in time. This is strictly equivalent to follow occupied states. 

Note finally the following useful property: for an independent particle state, = p. Indeed, 
the operator given by Eq. (jl7p then writes 

N 

P = Y1 (20) 

n=l 

Therefore, p is nothing but the projector on the subspace of occupied single-particle states of 
{(pui-ujsi) ■ Considering a set of orthonormal single-particle states, we have 

N 

= ^ Wm) WmWn) Wn\ = P- (21) 
m,n=l 

Assuming the above relation for the one-body density is equivalent to assume independent particle 
states. 

4- Dynamical Evolution from the Erhenfest theorem 

As discussed above, the basic approximation is to focus on a specific class of observables, which 
are the one-body observables in our case. Considering any observable F, using (Eq. ([2])) and its 
Hermitian conjugate for any N-body state {ip), we have 

= (V'l i?) F 1^) + (^1 F (^:^H =^{[H, F] )^ (22) 
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which is nothing but the standard Ehrenfest theorem for the operator F. The above equation is 
exact and vahd for any (correlated or not) state and observable. Here, H is the full microscopic 
Hamiltonian (see Eqs. @ and ([5])). For a one-body operator, F can be written as a linear 

combination of the a\aj. Therefore, it is sufficient to follow the expectation values of the {a\aj) , 
which are nothing but the matrix elements of p and are given by 

^ (alflj)^ = ^ Pji = «1 «i] )i>- (23) 

B. Time-Dependent Hartree-Fock (TDHF) 

1. Exact evolution of p^^"^ 



Our goal is to provide the best description as possible of the one-body density matrix evolution. 
The only approximation that will be made is that the system remains in an independent particle 
state at all time. Several derivations of TDHF exist so far. Here, we use the Ehrenfest theorem as 
a starting point. 

Reporting the Hamiltonian expression (Eq. dSJ) in the evolution of the pji (Eq. ([23l) ). we get 

^ ^ 0^ = ^kiiial aj , 4 ai] )v + ^ ^Mmn ( [aj aj , aj, a} a„ cim] (24) 

kl klmn 

Let US start with the first term associated to kinetic energy. Using Eqs. (IA4j) and ()A5p . we obtain 

{[alaj , alai])^ = {a\ aj alai)^ - {alai a] aj)^ 

= 6jk (aj ai)^ - (a| aj ai)^ - dn (dj, aj)^ + (dj. d| d; aj)^ 

= Sjk Pii - Pjk- (25) 

The kinetic energy term reduces to 

tki{[al aj , al ai] = {tjk Pki - tki Pjk) ■ (26) 
kl k 

For the two-body interaction, we should first express the expectation value of the commutator 

{[o-laj , alalanam])^, = {a] aj ala] anam)'^, - {ala] anamol aj)^ 

= {aj aj ttn cLm)^ Sjk - {al al an am)^, Sji + {al al aj aj an d^)^ 
-{al aj CLn CLj)^ 5nii + {al a] dm aj)^ Sni - {al a] al an dm d^)^. 

(27) 

The two terms with 6 annihilation/creation operators cancel out. Other terms are nothing but 
components of the two-body density matrix (defined in appendix IC]) . The two-body density can be 
decomposed into a sum of an anti-symmetric product of two one-body density (the uncorrelated 
part) plus the so-called two-body correlation matrix, denoted by C^^^ (see appendix[D]). Using this 
decomposition (Eq. (jPip ). Eq. ([27|) writes 



d| cij , al Oq cifi dm 



— {PmiPnl PmlPni ~l~ Cninil) ^jk ~l~ {PmkPni PmiPnk ~l~ Cnmik) ^jl 
+ {PjlPnk — PjkPnl + Cnjkl) 4ni + {PjkPml — PjlPmk + Cmjlk) 4i 

(28) 
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where exponents '^^^ and in equation (jPip have been omitted for simphcity. Altogether, the 
two-body interaction contribution to the one-body evolution reduces to 

4 ^ Vklmn a| CLj , aj. a| CLnCLm = - ^ [Vjklm {PliPmk — PlkPmi + Qmik) 

klmn klm 

+Vklim {PjlPmk — PjkPml + Cmjkl)] 



E 

klm 



Vjklm PliPmk + -^l mik 



-Vklim [ PjkPml + T^L^j mkl 



(29) 



where we have used Vkimn = —Vkinm = —vikmn and the anti-symmetry of C^^^ ^. The two-body 
contribution can finally be written as 

{U[p]jk Pki - U[p]ki Pjk) + - y^ {Vjlkm Ckmil — Vklim Cjmkl) (30) 
k klm 

where U[p] is the Hartree-Fock self-consistent mean- field. The latter can be written with the use 
of partial traces as 

U[p]ij = Yl ^^kji Pik = {i\ Tr2{v{l, 2) p{2)} \j) = Tt2{v{1, 2) p(2)},,. (31) 

kl 

The trace is only made on the second particle "2" while the labels i and j correspond to the particle 
labeled by "1". The matrix U[p] therefore corresponds to a one-body mean-field operator. In the 
following, we will often use the notation vi2 = v{l,2) and p2 = p(2). Then, Eq. (jSip simply reads 
U[p]i = Tr2{'Ui2 P2}- 

The second term in equation (jSOp reflects the effect of correlations on the evolution of one-body 
degrees of freedom. It can also be written with a partial trace on the particle "2" 

'Y^Viikm Ckmjl = Tr2{vi2 Ci2}ij- (32) 

klm 



Then, using Eqs. (^^, (p^ . (j26j) . (j30|) and (j32|) we finally deduce that the most general expression 
of the one-body density for any correlated system with a two-body interaction can be written as 

d 

"^^ -qI Pji = X] + Ujk[p]) Pki - {tki + Uki[p]) Pjk] 

+^ [ Tr2{-ui2 Ci2}ji - Tr2{Ci2 vi2}ji] ■ (33) 

The anti-symmetry of V12 and C12 implies that Tr2{!;i2 C12} = 2 Tr2{wi2 C12} and finally leads to 
the more compact form 

d 

^ ^ -Ql Pji = [h[p], P ]ji + Tr2 [ vi2, C12 ]ji , (34) 

where 

h[p]=t + U[p] (35) 

is the matrix associated to the one-body Hartree-Fock Hamiltonian. Note that, up to here, since no 
approximation has been made, the dynamical evolution of the one-body density is exact. However, 
the equation above requires a priori to also follow the two-body correlations in time which may be 
to complicated. 



The fact that C''^' is anti-symmetric is a consequence of the anti-symmetry of p'^' and can be deduced from 
anti-commutation rules for fermions (|A4|I entering in the definition of p'^' (Eq. ([C2}). 
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2. Time- dependent mean-field approximation 



In TDHF, two-body correlations are neglected at all time, i.e. C^^^ = O.This is equivalent to 
assume that the system remains in an independent particle state at all time (see appendix IF]) . 
Starting from Eq. (j34[) and neglecting the correlation term, we finally get the TDHF equation for 
the one-body density matrix p 



(36) 



In this section, we will now concentrate on the latter equation which could also be written in terms 
of operators as 



h[p\ , p 



(37) 



where h and p act both on the Hilbert space of single-particle states. In a complete basis of this 
space {|i)} with the closure relation \ = 1, we have 



p\i) = ^ Pji \ h\i) = h[p]ji \j) 



(38) 



Let us now give some properties of equation (j37p . First p = p is preserved. Therefore at all time, 



the density could be decomposed on a set of single-particle states p = |^'ra)(^'n| (Eq. (pO]) ). 

The TDHF equation can equivalently be written in terms of N coupled self-consistent equations 
on the single-particle states 



d 

ih g-^Wn{t)) = h[p{t)] Wn{t)), l<n<N. 



(39) 



Indeed, starting from the expression of p and using the above Schrodinger-like equation for single- 
particle states, we recover the TDHF equation (Eq. ([37]) ) that 



d ^ 
ih—p = ihy 
dV ^ 

n=l 



[kp\. 



This shows the equivalence of the single-particle representation and density formulation. 

Though Eqs (j39p take the form of Schrodinger equations, they are non llinear because the 
Hamiltonian depends on the one-body density. As a consequence it depends explicitly on time. 
We clearly see here some difficulties of mean-field theories. Indeed, we want to describe the system 
evolution between the initial and final time with a Hamiltonian which is itself depending on the 
evolution. 



C. The Skyrme effective interaction 



Mean-field equations have been derived for any general two-body Hamiltonian. However, in 
practice, the interaction is chosen to simplify numerical aspects. Most (i f not al l) applications of 
TDHF in the nuclear context have been performed using the Skyrme like [Skv56l | interaction. The 
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most widely used Skyrme force writes 



v{l,2) = to (1 + xoP^) 5 

1 
2 



+ ^ti [l + XlP^] ( k2 5 - 5 P 



+ iWo a - (kx 6kJ (40) 

where 5 = (5(f (1) - f (2)), k = (p(l) - p(2)) /?i (relative impulsion), a = <t(1) + <t(2), R = 
(f (1) + r(2)) /2, a{i) = dxii) e^+dij^i) ey+az{i) e^, (5"a^./y/^ (i) are operators acting on the spin of par- 
ticle i and are given in terms of Pauli matrices acting on the spin space. Pq- = (1 + <t(1) • <t(2)) /2 
corresponds to the exchange of the spin. p{r) = p^^\rsT, tst) is the particle density at r. "ti" 
and "t2" terms are non-local in space and simulate the short range part of the interaction. Finally 
"Wo" is the spin-orbit term. 

The very interesting aspect of this interaction is its zero range nature, which greatly simplifies 
the mean-field expression in coordinate space. Parameters {to-3, xo-3, Wo and a) are generally 
adjusted to reproduce nuclear properties like satura tion, incompres s ibility . .. of nuclear matter and 
selected properties of finite nuclei (see for instance Cha97l . ICha98l . lMevOCll | ) . 



An important aspect of the fitting procedure which has direct implication on nuclear reactions 
calculations is the following: fo r nuclea r structure calculation, center of mass contribution are 
removed (see section 3.2. of ref. Cha98l |) to obtain a better description of nuclei in their intrinsic 



frame. For collisions, only the intrinsic frame of the total system is considered and the same 
correction used for a single nucleus could not be used anymore. This is the reason why specific 
forces w hich explicitly do not account for center of mass correction have been developed like SLy4(i 
KimO 'il where parameters of the force have been adjusted to reproduce the same properties as 



Sly4 [Cha98] except that the center of mass correction is neglected. 

Note finally that the "is" term in Eq. (j40p explicitly depends on the system density. For 
specific integer values of a, the density dependent two-body interaction could be interpreted as 
a higher order interaction. In practice, non-integer values of a turns out to be more effective in 
reproducing nuclear properties. This however has important consequences. In particular, since 
the interaction depends on the system on which it is applied, strictly speaking, we cannot really 
use the terminology "interaction". One often use the very notion of Energy Density Functio nal 
to avoid confusion. In addition, translational invariance (but not Galilean invariance [Tho6|) is 
explicitly broken. 



D. Numerical implementation of TDHF: practical aspects 



tion 



Several a pplications of TDHF have been perforrned in th e last decade for nuclea r collec t ive mo^ 
1 studies f Sim03', NakOS', 'Alm05', Uma05, 'Rei06, ReiO TH and nuclear reacti ons fKimO?', SimOl 



Sim04 lUmaOGa . MarOG . Uma06b . Uma06c . Uma06d. .GuoOTl . lUmaOTl . lSim07 . SimOS . GuoOa 



Uma08l ] with various numerical methods to solve the TDHF equation. We present here a method 



used to implement TDHF for nuclear collisions. Different steps of a calculation are presented to 
better illustrate numerical constraints. To apply TDHF, we should 

• Construct the HF ground states of each of the collision partners. This requires to first solve 
their HF equations. 



• Starting from two Slaters at an initial distance Do, one should construct a single Slater 
associated to the composite system. 
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• The nuclei should be initially positioned and boosted to properly account for the reaction 
properties (Beam Energy, impact parameter...) 

• The dynamical evolution should be performed iteratively to solve the self-consistent TDHF 
equations. This is generally done by solving the Time Dependent Schrodinger equations on 
occupied states. 

• Finally, we should compute a set of observables to get informations on the reaction itself. 

Since all the information is contained in the one-body density, only one-body wave-function need 
to be considered. 

Sections III D 1 1 to III D 41 describe respectively the construction of initial Slater determinants and 
how initial conditions for reactions according to Rutherford trajectories are imposed. Section [IID 51 
presents the numerical implementation of mean-field transport equations. 



1. Hartree-Fock initial state 



We assume that the collision partners are initially in their ground states. Consistently with 
the TDHF approach, we should consider that the ground state of both nuclei are solution of the 
self-consistent Hartree-Fock (HF) states. The one-body density associated to each state is solution 
to the stationary version of the TDHF equation (Eq. (|37|) ): 

h[p] ,p\=0. (41) 

This equation is valid in any basis. Then, a specific basis should be chosen to explicitly solve the 
equation. Since h[p] and p do commute, we can choose common eigenstates. We denote this basis 
by {\a)} with h \a) = Cq, \a) and p \a) = Ua \a) ■ Here occupation numbers verifies tt-q, = ou 1. 
The A'^-body Slater determinant state is constructed from the occupied states (with = 1). 

Eigenvalues Cq, of h can eventually be interpreted as single-particle energies iVau72 l]. IntheHF 
approximation, the ground state is obtained by filling the lowest energy single-particle states. 
Therefore, we only need to find the N lowest eigenstates of h. 

2. Imaginary-time method 



The imaginary-time method Dav8Cll | is a widely used method to find the lowest eigenvalues of 



an operator (whose eigenstates is bound from below). Let us illustrate this method for a particle 
in an external field. The method consists in starting from an initial wave- function |fe) which is not 
a priori an eigenstate of the one-body Hamiltonian, denoted by H. This state can be decomposed 
onto the true eigenstates {H\pn) = En\pn)) 



^ ^ Cn\pr, 



(42) 



We apply the operator e on the initial state 



The lowest energy eigenstate \pq) associated to Eq can then be obtained from 
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Indeed, we have En > Eq and then Q-P{En-Eo) _^ g ^ _> qq except for n = 0, i.e only the 

ground state component does not vanish. The denominator is required because e~^^ is not unitary 

(/9 E M). The terminology "Imaginary-time" comes obvisouly from the fact that e~^^ looks like 

the propagator in time e~*^*/^ where the time is a complex quantity t = —ihf5. It is finally worth 
to mention that the initial state should be guessed in order to contain at least a small fraction of 
the lowest eigenstate. 

The imaginary-time method should be further improved to obtained eigenstates of the single- 
particle HF energies. Indeed, two additional difficulties exist: 

• Since the system is composed of N particles, we need the lowest eigenstates of h and 
therefore we need to guess initial starting points for the single-particle wave- functions. 

• The Hamiltonian h is non-linear and depends explicitly on the system density p. 

To solve the first difficulty, we apply the imaginary time method imposing that the states remain 
orthonormal (through a Graham- Schmidt orthogonalization for instance). We then expect to 
converge towards the lowest eigenstates if the states are not initially orthogonal to the "true" 
eigenstates. A common choice for the initial states are those of Harmonic or Nilsson potentials. 

To solve the second difficulty, we procede iteratively with small imaginary-time steps A/3. At 
each step, the density is calculated from the states and the mean-field Hamiltonian is modified 
accordingly. 

A schematic representation of the different steps of the imaginary-time procedure is 

{ki"V--|'^?)} P^"^ ^ = A/i[p(")] + (1 - A)/i(") 

^ ^ . (45) 

(n+l)v _ 1 /, ,v /„(n+l)|„/\ |„(«+1) ~ 



\V ') = irr \ W. 



7n 



Here A/i denotes the normalization of the state "i". Note that, "A" is an extra parameter generally 
used to slow down the convergence and avoid numerical instability BonOSl ] . In practice, it has 



the effect to mix the mean- field Hamiltonian at a step "n" with the one at the step "n-1". The 
small time-step increment are generally performed using simple development of the exponential, 
i.e. e-^^'*'"' ~ (1 - A/?/i(")). A gain, since the latter is not unitary a Graham- Schmidt method is 
used at each iteration. 

In practice, a specific Hilbert space basis should be first chosen to numerically express the 
different steps depicted in (j45p . The most common choice is either the space (resp. momentum) 
coordinate representation |rsr) ( resp. |pst)) or harmonic oscillator basis \nljm). In both cases, 
since an infinite number of states could not be considered, the basis has to be truncated. In 
coordinate case, this is achieved by restricting the space to a discrete finite box while in the 
Harmonic oscillator basis a limited number of major shells is considered. Once the basis is selected, 
all ingredients of the theory are expressed in this basis. 

At the end of the iterative procedure, a set of A'" single-particle states is obtained that completely 
defines the HF ground state. Unoccupied levels could also be obtained in a similar way although 
they do neither affect the density nor the mean-field. In Figure[Tl individual densities Yls l'/'i^('^)'^^P 
associated to neutron states in ^^O and obtained with the imaginary-time method are presented. In 
that specific case, the discretized coordinate space was used. Technically, the box should be taken 
sufficient large to avoid the effect of the boundary. Here Hard Boundary conditions are retained 
imposing that wave-functions cancel out outside the box. Note finally, that in Fig. [H only the 
lsl/2, lp3/2 and lpl/2 single-particle states are occupied ^. 



^ Each state presented here is doubly degenerated due to the exphcit assumption of time-reversal invariance in the 
calculation. 
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Figure 1: Spatial densities of single-particle neutron states in ^^O. The state lpl/2 is cut at the middle to 
show the hole in its center. 



3. The two nuclei case: initial state construction 

Since each nucleus has been constructed in its HF ground state, we start from two independent 
particles systems to construct the initial TDHF condition. However, we have shown previously 
that Eq. (j36p describes a priori the evolution of a single Slater determinant. 

It is however always possible to construct a Slater determinant = associated to 

A = Ai + A2 independent particles from two Slater determinants = l^a^-.-o^^ ) and |02) = 
\4>Pi---I3a2 ) systems are initially well separated. Let us consider the two one-body densities 

pi = X^nii |Q^n)(an| and p2 = 'Y^n=i For each densities we have pf = pi (see section lTl A 3p . 

If we now construct the square of the total density p = pi + P2 

= Pi + P2 + h h + h Pi= P+^^^^ (|am)(am|/?n)(/3n| + H.c.) (46) 

m=l n=l 

where H.c. stands for Hermitian conjugated. Therefore, to have fP = p, single-particle states should 
not overlap. In practice, this is possible due to "Hard boundary" conditions in the HF calculations 
which impose that the single-particle wave-functions vanish at some distance. Therefore, we should 
just be careful to consider a box for TDHF that is large enough to avoid overlap between the two 
boxes used in the HF case (see figure [2]) . 



4-. Dynamical evolution of nuclei 

The TDHF theory is a quantum theory since it treats explicitly the particles with wave- 
function. However, the restriction to independent particles states does not allow, in general, for a 
probabilistic interpretation of reactions channels. Accordingly, some quantum aspects are missing, 
we will clearly see this pathology in fusion reactions (where the fusion probability will either be 
or 1). Essentially, TDHF gives classical trajectories for the evolution of centers of mass. 

We consider the system in its total center of mass frame. The impact parameter and beam 
energy (which fixes the initial velocity v of the projectile at infinite distance) allow us to determine 
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Figure 2: Schematic representation of the 2 initial HF boxes positioned in a bigger box for TDHF calculations. 



the initial TDHF condition. At this initial time (t = 0), the two nuclei have a relative distance Dq. 
In general, a Rutherford trajectory is assumed to account for Coulomb trajectory from infinity to 
Do- This hypothesis is consistent with the assumption that the two nuclei are in their ground state 
at t = 0, i.e. we assume that no energy has been transfered from the relative motion to internal 
degrees of freedom up to Dq. 

Using notations of section HID 31 a velocity Vj is applied to each nucleus (i = 1 or 2) imposing 
the im pulsion Pj = Ajmvj. This is performed, by applying a translation of each pi in momentum 
space |Tho62l | 

= 0) = e*™^'-*^/^ Pi e-'""^^-^/^ (47) 

where the position operator r = x ex + y Gy + z ez acts in single-particle space. Note that, here pi 
is reserved to the static HF while pi{t = 0) with time in parenthesis corresponds to boosted HF. 

In practice, since one usually follows directly single-particle wave- functions, the translation in 
momentum space is directly applied to the waves 

^^Jr,t = 0) = e^— ^ ifl^ir) 1 < n < 

^l{r,t = 0) =e*™v2.f ^.tjr) l<n<A2. (48) 

Once the two nuclei are positioned on the network and properly boosted, there is no more reason 
to distinguish single-particle states from one or the other collision partner. 



5. Numerical methods for dynamics 

To solve the system of equations of motion, we should solve TDHF equations for occupied states 
(Eq. (j39|) ). The main difficulty is the fact that the Hamiltonian itself depends on time. As a 
consequence, as in the imaginary-time case, specific procedure should be implemented to account 
for self-consistency in propagators. We consider a small time step increment At and perform 
the time evolution iteratively. Over small time intervals At, the Hamiltonian is almost constant. 
However, to conserve the total energy, the numerical algorithm should be symmetric with respect 
to time-reversal operation. This implies to consider the Hamiltonian value at time t + ^ for the 
evolution of wave-functions from t to t + At Flo78l | ^ 



\u{t + At)) ^e-'^~^i'+^) \u{t)). (49) 



^ This algorithm is similar to a Runge-Kutta method. 
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A schematic illustration of the real time propagation could be written as: 



{|-i"V->?)} 

il 

(n+l)\ -iMh("+i'> I (n)v 
A >)=e 'n"" \ul ') 

it 



.(n) 



— 2 



p{n+l) 



(50) 



where corresponds to an approximation of \^{tn = nAt)). In this algorithm, starting from 

the density at time t, a first estimate of the density at time t + At, denoted by obtained. 
The Hamiltonian used in the propagator is then computed using the average density obtained from 
p^") and p^"'~^^\ Then, the real new density at time i + At is obtained using this Hamiltonian. 
As in the imaginary-time case, an approximate form of the exponential is generally used which in 
some cases, breaks the unitarity (even in the real time evolution) and orthonormalization of the 
single particle states must be controlled. 



III. APPLICATION OF TDHF TO REACTIONS AROUND THE FUSION BARRIER 



TDHF has been applied to nuclear ph ysics na o re than thirty years ago. First applications were 
essentially dedica ted to fusion reactions Bon76l . lBon78l . lFlo78l . lNeg82i] . Figure El adapted from 



P. Bonche et al. Bon78l ] . illustrates the predicting power of TDHF for fusion cross sections. The 



main difference with actual calculations is that, at that time, several symmetries were used to 
make the calculation tractable. The major advantage of imposing symmetry was to reduce the 
dimensionality (and therefore the numerical effort). The major drawback was the reduction of 
applicability. The second difference with nowadays calculations comes from the fact that simplified 
forces were used, missing the richness of effective interactions used in up to date nuclear structure 
studies. Recently, all symmetry assumptio ns have b een relaxed and forces using all terms of the 
energy functional have been impleme nted Uma06al ] . We present here several applications in 3D 
coordinate space with the full SLy4(i fKimOTf Skyrme force. 

First, general aspects related to fusion reactions are presented, like fusion barrier properties, 
cross sections... Then, we illustrate how a microscopic dynamical model can be used to get infor- 
mations from the underlying physical process. Besides fusion probability, transfer of nucleon will 
be discussed. Finally, limitations of standard mean-field models are discussed. 



A. Selected aspects of fusion 



1. Definition of fusion 



Nuclear fusion is a physical process where two initially well separated nuclei collide and form a 
compound nucleus which has essentially lost the memory on entrance channel (Bohr hypothesis). 
This hypothesis is rather simple while the experimental measurement is rarely trivial. Indeed, 
energy and angular momentum conservations imply compound nucleus generally formed at rather 
high internal excitation and angular momentum. As a consequence, the fused system cools down by 
gamma, particle emissions and/or eventually fission. Overall, several decay channels are competing 
leading to a broad range of final phase-space which eventually overlaps with direct reactions (and 
more generally pre-equilibrium) processes. This is for instance the case of quasi-fission (where the 
systems keep partial memory of the entrance channel) which leads to mass and charge distribution 
which can eventually be similar to fission. Similarly complete and incomplete fusion are sometimes 
difficult to disentangle due to the presence of direct break-up channels (which are enhanced in the 
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Figure 3: The ^^O+^^O fusion was one of the first application of TDHF calculation (adapted from Bon78| ) 



weakly bound nucleus case). Therefore, fusion events are sometimes experimentally difficult to 
distinguish from other processes. 



2. One dimensional approximation 

The simplest approach to fusion reaction is to consider the relative distance r between the 
centers of mass of the nuclei as the most relevant degree of freedom. Fusion reactions then reduce 
to the dynamical evolution in a one-dimensional potential where the potential is deduced from 
the long-range Coulomb repulsive interaction of the two nuclei and from their short range mutual 
nuclear attraction (see figure |4]). We assume that fusion takes place when the system reaches the 
inner part of the fusion barrier (r < rjy on figure H]). Within this approximation, fusion cross 
sections can be expressed as a function of the transmission probability Ti{E) for each energy E 
and angular momentum I 



ajUE) = (2£ + l)r,(i?), 



(51) 



where denotes the reduced mass. Below the fusion barrier, i.e. {E < i?), fusion is possible only 
by quantum tunneling. Using the WKB (Wentzel-Kramers-Brillouin) approximation leads to the 
following transmission coefficients 



1 -|- exp 



dryj2ii{U{i,r) - E) 



rc 



(52) 



where f7(£, r) is the total (nuclear -|- Coulomb -|- centrifugal) potentials while rc and correspond 
to "turning points" at energies E (see Fig. H|. Another simplification can eventually be made using 
a parabolic approximation. In that case, the potential is approximated by an inverted parabola 
with curvature Hlob- This approximation is justified close to the fusion barrier only and leads to 
analytical expression for the transmission coefficients 



Ti{E) 



1-hexp [2t:{B - E)/nu:B] 
Finally, summation on different I leads to the Wong formula jWon73| | 



crfus{E) 



fllOB 



2E 



In < 1 -|- exp 



27r 



(E-B) 



(53) 



(54) 
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Figure 4: Left: example of nuclear, Coulomb, centrifugal (with I — 40), and tot al = and ^ — 40) 
potentials for ^®0+^°^Pb. The Wong parameterization is used for the nuclear part [Won73j . Right: zoom 
around the barrier for i — 40. For a given energy i?, the inner and outer turning points at r^v and rc 
respectively, as well as the parabolic approximation (dashed line) are presented. 



Predictions of this one-dimensional approximation are compared to experimental observations for 
0+Sm and Ni+Ni systems in Fig. [5j The comparison is relatively satisfactory above the barrier 
but fails to reproduce sub-barrier cross-sections which are underestimated by several orders of 
magnitude by the Wong formula. 

In addition, experimental data clearly show large differences from one isotope to the other which 
underlines the inherent nuclear structure effects and could not be simply explained by the change 
of the radii. Note finally that improvements where the parabola and/or the WKB approximations 
are not made do not improve significantly the comparison. 



3. Coupling between relative motion and internal degrees of freedom 

The discrepancy between experiments and simple estimates can directly be traced back to the 
fact that we assumed nuclei as rigid objects without internal structure. In fact, fusion is affected by 
the reorganization of internal degrees freedom as the two nuclei approach. This induces a coupling 
between the internal degrees of freedom and the relative motion. The only way to include this effect 
is to add additional degrees of freedom in the description of fusion. At the macroscopic level, this 
is generally achieved by introducing, for instance, deformation, orientation, neck parameters, mass 
and charge asymmetries leading to more complex "multi-dimensional" potentials (see Lac02[ | and 



references therein). This will be illustrated with TDHF calculations below, in particular to study 
the effects of deformation. 

Sub-barrier fusion is a perfect illustration of effects induced by couplings to internal degrees of 
freedom. These couplings can lead to an ensemble of fusion barriers called harriers distribution 
and can give enhancement of cross sections by several order of magnitude in the sub-barrier fu- 
sion regime^. A proper description of these effects could only be achieved if inelastic excitations. 



Note that, as a counterpart, fusion above the barrier are usually reduced compared to the one-dimensional case. 
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Figure 5: Fusion cross sections for different systems as a function of center of mass energy (divided by their 
respective fusion barrier energy). Left: ^^O projectile on different Samarium isotopes. Right: Nickel on 
Nickel collisions. Curves correspond to calculations deduced fro m barr ier penetration model with standard 
parameterization of nucleus- nucleus potentials (extracted from Das98l |). 



collective modes, transfer and all relevant processes are properly accounted for (see Bal98l . Das98l ]). 

Note finally that additional important effects could appear when weakly bound nuclei are in- 
volved. In particular, break-up channels and new collective modes may become important. Prom 
an experimental point of view, influence of these new effects on reduction/increase of fusion cross 
sections is still under debate. With future low energy radioactive beams, we do expect to get 
additional informations on the reaction mechanisms with weakly bound nuclei around the barrier. 



Barriers distribution 



The experiment al barri ers distribution D'^^{E) is obtained from the excitation function afusiE) 



using the relation [RowQll ] 

dE2 



D^iE) = (55) 



This function could be interpreted as the probability that the system has its barrier at energy E. 

Let us illustrate this with a simple cases. Consider a classical model with a single barrier at 
energy B. In that case, Db identifies with the Dirac function 6{E — B). The latter formula is 
consistent with the fusion probability found in classical systems (i.e. /i ^ 0) deduced from the 
Wong formula (j54p . The fusion probability is zero below the barrier and equal to 7ri?^(l — ^) for 
E > B. We therefore deduce that the distribution barrier reads Db{E) = 7tR^6{E — B) which is 
the expected behavior. For quantum systems, Db{E) spreads over a wider range of energies. 



B. Fiision barriers and excitation functions with TDHF 



In this section, TDHF results essentially extracted from Refs. 



SimOSl . IUma06b | are presented 
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Figure 6: Schematic illustration of separation plane between two fragments after the contact. 



1. Fragments trajectories 

A way to characterize TDHF trajectories consists in matching the microscopic theory with the 
one-dimensional model described previously. The main difficulty is then to define the relative 
distance between the two nuclei. When fragments are well separated, such a distance could be 
easily defined. On the opposite, after the touching, this becomes more complicated and could only 
be achieved for short time after the touching. Fortunately, around the barrier, nuclei are generally 
still well separated. In practice, the main ax is of th e reaction could be extracted (it generally 
corresponds to the main axis of deformation) jWasOSl ] . Then, along this axis, the local minimum 
of the density profile defines the separation plane (see figure [6]) . Note that this distance does not 
cancel out even for one sphere. 

2. Barriers distribution Jar two initially spherical nuclei 

We first consider ^^O+^'^^Pb reaction. The two nuclei are initially in their ground state and 
are at a relative distance Dq = 44.8 fm. Evolutions of relative distances for head-on collisions 
(6 = fm) obtained from different center of mass energies (from 74.2 et 75 MeV) are displayed in 
Fig. [7|^. Due to the narrow range of initial center of mass energies, relative distance evolutions 
before touching are all similar. On opposite, after t = 600 fm/c, we clearly see two groups of 
trajectories: 

• E < 74.44 MeV, the two fragments re-separate. 

• E > 74.45 MeV, the relative distances remain small (r < 10 fm). 

The latter case corresponds to the formation of a compound nucleus after the collision. A precise 
analysis shows that the predicted fusion barrier lies between 74.44 and 74.45 MeV. Experimentally, 
it is found around 74 MeV, with a width of 4 MeV. It is rather interesting again to mention that 
TDHF seems to precisely describe fusion barrier while no parameters of the effective interaction has 
been adjusted on reactions. In order to better understand phenomena occuring in TDHF around 
the fusion barrier, the local part of the one-body density (p(r) = X^^,- /o(rsr, rsr)) evolutions are 
displayed in Figs. [8] and [9] respectively for energies just below and above the estimated fusion 



Each trajectory presented in figure [7] have been computed using 4 hour CPU time on a NEC/SX-8 processor. 



25 




Fusion 

I I I I I I I I I I I I I i_ 

200 400 600 800 1000 1200 

time (fm/c) 

Figure 7: Relative distance between fragments as a function of time for head-on ^^0+^°*Pb reactions. 




Figure 8: Density evolution for the reaction ^^0+^°*Pb corresponding to head-on collision at center of mass 
energy 74.44 MeV (just below fusion barrier). The red surfaces correspond to an iso-density half of the 
saturation density. Each figure is separated by a time step of 135 fm/c. 



barrier. In tiie former case, the composite system hesitates to fusion. It forms a "di-nuclear" 
system for relatively long time (~ 500 fm/c) before re-separating. During this time, nucleons are 
exchanged between the two partners. In the second case, the system passes the fusion barrier (it is 
just 10 keV above the barrier). More generally, the two figures illustrates the richness of physical 
phenomena contained in TDHF: surface diffusivity, neck formation at early stage of fusion process, 
quadrupole/octupole shapes of compound nucleus... 

A similar agreement between experimental and calculated fusion barriers is found in other 
systems as shown in Fig. [TUJ Several projectiles starting from light ^^O to medium mass ^^Ni 
nuclei and targets from ^''Ca to ^38^ have been considered. The lowest energy barrier corresponds 
to "^"Ca-F^^Ca while the highest is obtained for ^^Ti-F^ospb. We clearly see in this figure that 
barriers extracted from TDHF, where the only inputs are the ef fective forces p arameters |Kim97l | . 
give a better agreement with data than Bass empirical barriers |Bas77l . feasSOl ]. 




Figure 9: Same as figure [8] for center of mass energy 74.45 MeV (just above fusion barrier). 
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Bevp (MeV) 

Figure 10: Macroscopic fusion barriers from the Bass parametrization (stars) compared to barriers extracted 
from TDHF (circles) as a function of experimental barriers (centroids of fusion barrier distributions except 
for the last point corresponding to ^'*Ti+^°*Pb_where the barrier distribution is obtained from quasi-elastic 
scattering (MitOTj l. Fi gure extracted from [Sim08l |. 
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Figure 11: Typical experimental barrier distributions obtained for a light spherical nucleus on a heavy 
deformed one. Arrows indicate barriers deduced from TDHF assuming different orientations of the deformed 
nucleus at the touching point. 



3. Barrier distribution from collisions between a spherical and a deformed nucleus 

Fusion involving at least one deformed nucleus is helpful to illustrate the appearance of several 
fusion barriers which could be easily interpreted in terms of a classical variable, the "relative 
orientation" of the two nuclei. Figure [TT] gives examples of experimental barriers distribution 
(Eq. (j55p ) for reactions involving one spherical light nucleus and a heavy deformed one. These 
distributions have a typical width of ~ 10 MeV. Tunneling effect can generally account for ~ 2 — 3 



MeV widths Row9lf | . An alternative interpretation should then be invoked to understand such a 
spreading. 

One can understand this spreading as an effect of different orientations of the deformed nucleus 
at the touching point. Indeed, for an elongated nucleus, if the deformation axis matches the reaction 
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axis (at 6 = fm for instance) , then the barrier is reduced due to the plug in of nuclear effects at a 
larger distance. On opposite, when the deformation axis is perpendicular to the collision axis, the 
barrier is increased. This two extreme cases correspond to an upper and lower limits (indicated by 
arrows on Fig. [TTI) for the barrier distribution. 

This effect is illustrated in figure [12] where potentials extracted from TDHF calculations are 
displayed as a function of relative distance for two different initial orientations of a prolate nucleus. 
In this case, when the orientation and collision axis are parallel, the contact between nuclei takes 
place at larger relative distance as compared to the perpendicular case which leads to a more 
compact configuration. The fusion barrier is then reduced for the parallel case because the Coulomb 
repulsion is smaller when the nuclear attraction starts to be significant. We also note in Fig. [11] 
that barriers distributions are peaked for compact configurations. This means that the probability 
that the system goes through a high energy barrier is higher than the one associated to low energy 
barriers. This is due to the fact that the deformed nucleus is prolate and not oblate. Let us consider 
for simplicity that the collision axis is the z axis and that the deformation axis could only take 
three orientations along x, y oi z with equal probabilities. For a prolate nucleus (the elongation 
occurs along the deformation axis), the compact configuration (associated to higher barriers) is 
reached if the deformation axis is a; or y, then with a probability |. Having a lower barrier is then 
less probable, corresponding to a deformation along z with a probability ^. The opposite occurs 
with an oblate nucleus (the elongation is perpendicular to the deformation axis). In this case, only 
an orientation along the z axis leads to a compact configuration, then with a probability i. This 
simplified discussion gives a qualitative understanding of figure [11] 

In the simple discussion above, we have assumed that the distribution of relative orientations 
is isotropic. This hypothesis generally breaks down due to the long range Coulomb interaction 
which tends to polarize the nuclei during the approaching phase. In the case of a prolate nucleus, 
the Coulomb repulsion being stronger on the closest tip of the deformed nucleus, the net effect 
is to favor orientations where the deformation axis is perpendicular to the collision axis. This 
polarization effect modifies the simple isotrop ic picture, in particular when the deformed nucleus 
is light and its collision partner heavy SimOj]. 

In summary, we have seen that barrier position and height can be affected by the structure 
of the collision partners. It is hazardous to conclude that fusion could be a tool to infer nuclear 
structure properties. However, with precise fusion measurement, one can clearly get informations 
on deformation properties which might be hardly reached with other classical techniques of spec- 
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Figure 13: Experimental fusion cross section (stars) compared to cross section deduced from TDHF (lines) 
for ^^0+^"'*Pb. The two lines correspond respectively to lower and upper limits of theoretical cross sections. 



troscopy. 



4- Excitation functions 

The calculations presented in the previous section were performed using TDHF at zero impact 
parameter. To compute excitation functions (Fig. [3|), calculations should include all impact pa- 
rameters up to grazing. The fusion cross section is given by Eq. (j51|) where T(^{E) is nothing but 
the fusion probability for a given center of mass energy E and angular momentum \Jl{l + l)fi. The 
independent particle hypothesis implies that T^i^E) = 1 for / < Imax. 

(E)and for / > lmax{E). 



This finally leads to the so-called "quantum sharp cut-off formula" |Bla54l ] 



^fus{E) = ^{lma.{E) + lf. (56) 

To avoid discontinuities due to the cut-off and integer values of ImaxiE), {lmax{E) + l)h is generally 
approximated by its semi-classical equivalent Cc = ^2^E be- The latter corresponds to the classical 
angular momentum thres hold fo r fusion and he denotes the maximum impact parameter below 
which fusion takes place BasSOl ] . This replacement is justified by the fact that {Imax 

+ Vf and 



Cl/h^ are both greater than lmax{lmax + 1) and lower than {Imax + l)(^max + 2). Accordingly, we 
finally obtain the standard classical expression for fusion cross sections (yfus{E) — tt£.'^/2iiE = irb'^. 

Figure [13] presents a comparison between calculated and experimental cross sections for the 
i6Q_j_208p]^ system. Discretizing the impact parameters gives an upper and lower bound for fusion 
cross sections. We essentially see that above the fusion barrier, TDHF calculations reproduce 
rather well the experimental observations (the cross section is however slightly overestimated by 
16%) while at energies below the Coulomb barrier, the calculation misses the quantum tunneling 
contribution. 
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C. Nucleon transfer around the fusion barrier 



In this section, we study nucleon transfer below the fusion barrier with TDHF. The basic 
observable associated to transfer is simply the mass (or nucleon number) of the two fragments 
after re-separation. If the latter differs from the entrance channel, this is an obvious signature of 
transfer [UmaOSi] . Another signature would be the variance of nucleon number in the fragments. 



1. Transfer identification 

As illustrated in figure El two nuclei can form a di-nuclear system with a neck and then re- 
separate. There is a priori no reason that these two fragments conserve the same neutron and 
proton numbers as in the entrance channel (except for symmetric reactions). Indeed, between 
the touching and re-separation, nucleons can be exchanged. In TDHF calculation, this exchange is 
treated through the time-dependent distortion of single-particle wave-function which can eventually 
be partially transfered from one partner to the other. 

The following operator written in r-space defines the number of particles in the right side of the 
separation plane (defined arbitrarily as x > 0): 

Nd = Y^ [dr d\r st) a(rsr) W(x) (57) 

ST 

where TC{x) is the Heavyside function equal to 1 if a; > and elsewhere. 

Denoting by {i\j)D = Z^sr/'^'' V'j^i^) "^{x) the overlap (limited to the right side) be- 

tween two single-particle states and using Eq. (jFSP for an independent particles state \4>), we obtain 
(in the specific basis which diagonalizes the one-body density associated to \4>)) 

{Nd)4, = {i\j)D (4 dj)^ = {i\i)D rii. (58) 

ij i 

Figure O gives the average final neutron and proton numbers of the smallest fragment in exit 
channels of ^^0-l-^'^'^Pb reaction as a function of center of mass energy. We see that the more the 
energy increases, the more the ^^O looses protons. At an energy just below the barrier, it has 
transfered around 2 protons to the ^'^^Pb. A possible explanation is the fast N/Z equilibration 
which is expected to take place at contact. Indeed, considering the case where exactly two protons 
have been transfered leads to the chemical equation 

1^0(1) +208pb(1.54) ^1^0(1.33) +2i0po(i.5) 

where quantities in parenthesis correspond to the N/Z values. We see that the N/Z initial asym- 
metry is almost equilibrated in the exit channel. This equilibration process also takes place 
at the f irst ins t ant of fusion, and is increased when N/Z differences between partners increase 
BonSll . ICho93l . ISimOll . ISimOTi ]. 



Note finally that these results on transfer agree qualitatively with what is observed experi- 
mentally Vul86l | . Indeed, it is observed that the one proton transfer channel from ^^O to ^'^^Pb 



dominates below the Coulomb barrier, while at higher energies the two protons transfer becomes 
the main transfer channel. 



2. Many-Body states associated to each fragment 



Well below the fusion barrier, transfer is prohibited and exiting fragments essentially refiect 
mass and charge properties of the entrance channel. In this case, the ^^O remains an eigenstate 
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Figure 14: Average neutron (solid line) and proton (dotted line) numbers of the lightest fragment in exit 
channel of head-on ^^0+^*'^Pb collisions below the Coulomb barrier as a function of the center of mass 
energy. 



of the proton and neutron number operators . For these state, the variance of particle number 

operator (in the right side) is strictly zero: ctd = \J {N'^) — {Nd)'^ = 0. 

This property is lost at higher energies where transfer occurs. Then each fragment in the exit 
channel does not have a "good" particle number and can therefore not be associated anymore to 
a single Slater determinant. It corresponds to a (eventually complicated) correlated state even if 
the wave- function associated to the total system is a Slater determinant. 

Le t us calculate the variance of Nd after the reaction. We follow here the technique of Dasso 
et al. [Das 79]. Using anti-commutation relations ()A5p . closure relation = i, and Eq. (jFSp . 

we get in the same way as for Eq. (j58p 

(Nl)^ = ii\j)D {k\l)D {al^j^ai)^ (59) 

ijkl 
ijkl 

= Y {Ai)D + Y ((^IOd 01j)d - K«Ij)d|^) ■ (60) 

The variance must then follow 

TV N 
i=l *)i=l 

Applying this formula to the small fragment in the exit channel at a center of mass energy 74.44 
MeV, we get a Dp — 0.5 for protons and a^^ — 0.3 for neutrons. This deviation from zero clearly 
indicates that the Many-Body state on each side of the separation plane is not a pure state with a 
good particle number anymore but a complicated mixing of states with various particle numbers. 
In addition, even if the neutron number is almost constant in average, this non-zero variance is a 
signature of neutron transfer. 



We may introduce the terminology that a Many-Body state has a "good" neutron (or proton) number when the 
state is an eigenstate of the particle number A'^ = 
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Note finally, that the variance is maximum when the overlap (restricted to the right side) 
between two wave-functions is zero, i.e. if {i\j)D oc 5ij. Then Eq. (j6ip identifies with a binomial 
distribution 



N 



1=1 

This gives an upper limit for the variance 

CTD < 



This upper limit is an intrinsic limitation of independent particle systems which will not be able 
to reproduce distribution with fluctuations greater than this limit. This is a clear limitation of 
TDHF, which in general un deresti r nates w idths of fragments mass and charge distribution observed 
in deep inelastic collisions [Das79l . [Goe82t | . The account for correlation beyond the single-particle 



approximation is the only way to escape this difficulty. 

D. Summary: success and limitations of TDHF 

In this section, we have compared TDHF calculation with experiments on the different following 
physical effects and quantities 

• Barriers positions and height 

• Shape of barrier distribution for deformed nuclei 

• Excitation function 

• Populations of fragments mass and charge below the fusion barrier. 

Overall, TDHF gives a reasonable qualitative (and sometimes quantitative) agreement with exper- 
iments. We have also pointed out some of the limitations of TDHF: 

• Essentially tunneling below the Coulomb barrier is not accounted for 

• Fluctuations of one-body observables like particle numbers are underestimated. 

Beside nuclear reactions, differences between TDHF calculations and experiments could also be 
observed in giant collective vibrations. For i nstance , damping of collective motion could not be 
properly accounted for at the mean- field level [Lac04i |. This could again be traced back to missing 



two-body correlations effects. The rest of this lecture is devoted to transport theories going beyond 
mean-field. 

IV. DYNAMICAL THEORIES BEYOND MEAN-FIELD 

In the previous sections, we have illustrated some successes of mean-field theories in the descrip- 
tion of nuclear dynamics. The application of TDHF to nuclear excitations or reactions was a major 
advance. At the same time, the mean-field theory neglects physical effects which play a major role 
in nuclei. For instance short and long range correlations in static nucl ei coul d only be accounted 
for by a proper treatment of pairing effects and configuration mixing Ben03t |. Conjointly, as the 



collision energy between two nuclei increases, the Pauli principle becomes less effective to block 
direct nucleon-nucleon collisions. Then, two-body correlations should be explicitly accounted for. 
The description of such two-body effects is crucial since, for instance, this would be the only way 
to understand how a nucleus could thermalize. During the past decades, several approaches have 
been developed to introduce correlations beyond mean-field. We summarize here some of these 
extensions. 
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A. Time-dependent mean-field with pairing correlations 



Pairing correlations are sometimes needed to des cribe nuclear systems. However, they are 
neglected in the independent particles approximation RinSOl . iBriOSl ] . Mean-field approaches could 
naturally incorporate this effect by considering more general Many-Body states as products of 
independent quasi-particles. This leads to the so called Time-Dependent Hartree-Fock-Bogoliubov 
(TDHFB) theory. 



1. Quasi-particle vacuum 



Wick's theorem (see appendix [A]) simplifies the calculation of observables expectation values on 
vacuas. In appendix|El we show, on one hand, that an independent particles state |0) is a vacuum 
(called HF vacuum) for a specific set of operators, and, on the other hand, is a simple example 
of "quasi-particle" vacuum. The canonical and unitary transformation connecting quasi-particle 
creation and annihilation operators (denoted by Pa an d 0„ re spectively) to particle ones is given 



by Eq. ()E4p . It is called Bogolyubov transformation |Bog58j]. The vacuum associated to these 
quasi-particle operators is given by 

k 

which obviously insures Pi\ipbogo) = 0. This state is then called a quasiparticle vacuum. 

Let us illustrate how pairing could be included in such a Many-Body state. We consider the 
specific case where the transformation is simply written as 

Pp ~ ^p^p ^p^^i Pp ~ ^p^p ~t~ ^p^p^ (^^) 

In this case, Vp corresponds to the occupation probability of the state p while Up = l — Vp. Eq. (j65|) 
is nothing but the standard BCS transformation jBar57t | where single-particle states p and p form 
a Cooper pair. We use the standard convention p = —p < 0. Starting from Eq. (|64p . grouping 
quasi-particles operators by pairs and using anti-commutation rules (|A4tiA5p . we can rewrite the 
Many-Body state as 

p>0 

The above expression corresponds to the standard BCS state expressed in terms of pairs of corre- 
lated nucleons {p,p}. In fact, any general Bogolyubov tra nsformation can be wri tten in a BCS form 
using the Bloch-Messiah-Zumino decomposition theorem Blo62l . [Zum62l . [Rin8Clf | . We therefore see 



that to account for more general Many-Body states is a way to incorporate correlations beyond 
mean-field. 



2. Expectation values of observables on quasi-particle vacua 

Similarly to HF vacua, the Wick's theorem applies to Quasi-particle vacua. In this case, however, 
contractions associated to the one-body density p^/B = {a^f^CLa) are not the only one which differ from 

zero. The anomalous density, with matrix elements Ka/B = {apcia) (which also implies k*^ = (aaajj)) 
should also be taken into account. The latter contractions cancel out for independent particle 
systems. The Bogolyubov transformation (Eq. ()E4p ) can be inverted to express the a) and a in 
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terms of quasi-particles operators {P\(3) : 



Using these expressions in p et k and the fact that only {iSi^j) differ from zero, we deduce 



Pa/S 



y*jjl 



(67) 



(68) 



These contractions are used to write the generahzed density matrix defined as 



n 



\[{d]4)) ({aA)) J 



p K 
-Hi* 1-1 



(69) 



The anomalous density enables us to treat correlations that were neglected at the independent 
particles level. Indeed, elements of the associated two-body correlation now read 



p\fki = {ij\pi2\kl) = {dlajaja 



t A. . ci^dj a|dj + d|,d| ajdi 



PikPjl - PilPjk + KijKli- 



(70) 



Countrary to Slater determinant, C12 does not a priori cancels out. We further see that the HFB 
theory leads to a separable form of the two body correlation 



ijkl — ^ij^kl- 



(71) 



In turns, HFB is more complex than the HF case. For instance, the state is not anymore an eigen- 
state of the particle number operator. The symmetry associated to particle number conservation 
is explicitly broken. Fluctuations associated to the particle number N = doda now write 



(iV2) - {Nf = 2 TY(KKt) = 2Tt{p- 



(72) 



In general these quantity is non-zero for a quasi-particle vacuum. The fact that the particle number 
is not conserved implies that it has to be constrained in average in nuclear structure studies (this 
is generally done by adding a specific Lagrange multiplier to the variational principle). It is worth 
mentioning that in a TDHFB evolution using a (possibly effective) two-body interaction (Eq. dS])), 

the expectation value of iV is a constant of motion. Therefore, no constraint on particle number is 
necessary in the dynamical case. 



TDHFB Equations 



D ifferent techniques cou ld be used to derive the equation of motion in the HFB approxima- 
tion [RinSd . lBen03l . pa86]. Here, the same strategy as section IIIBI is followed. The evolution 
of the Many-Body HFB state is given by the evolution of its normal and anomalous densities 
Pij = (ajdj) and Kij = (djdj) which are given by the Ehrenfest theorem 



in—K 
dt 



dV ^ 



d]d j,H 



H 



(73) 
(74) 
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a. One-body density evolution 
Eq. (j73p can be identified to Eq. (j34p which has been obtained for any correlated state and any 
two-body Hamiltonian. The only difference comes from the simplified correlations for quasi-particle 
vacua (Eq. (fTTI) ) 



klm 



(hjkPki — hikPjk — ^jkK-li + KjfcA^J 



(75) 



where we have introduced the pairing field 

Ajj = - ^ VijklKkl- 



kl 



(76) 



b. Evolution of k 
Similarly, the equation of motion for the anomalous density k reads 



d 



aiaj,d\ai 



1 



+ T ^ Vklmni 



kl 



klrnn 



— ^ '^ jtjk'^ki tik^kj) + ^ '^ {'^kjlmPlkf^mi VkUmPlkf^mj , 



+ ^ y^^jv kjlmPik^l 



klm 



m ^kilmPjk^lm , 



klm 



{hjk^ki + l^jkh*ki — ^jkPki ~ Pjk^ki) + Ajj. 



(77) 



Eqs. ([75]) and (j77j) give the evolution of the matricies p and k 

ih^p = + kA* - Ak*, 

dt 

ih—K = hK + Kh* - pA- Ap* + A. 
dt 



(78) 
(79) 



Finally, using the generalized density matrix TZ and the generalized HFB Hamiltonian 7i defined 
as 



n 



h A 

-A* -h* 



(80) 



Eqs. (ITSp and ()79p can be written in a more compact form 

at 



(81) 



The above TDHFB equation generalizes the TDHF case (Eq. (j36p ) by accounting for pairing effects 
on dynamics. 
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4- Application of TDHFB theory 

In practice, numerical implementation of TDHFB is much more complex than the TDHF case. 
This is certainly the reason why, although first applications of TDHF started more than 30 years 
ago, very few attempts to apply TDHFB exist so far. This could be traced back to conceptual 
and practical difficulties inherent to this theory. The first difficulty comes from the fact that 
TDHFB equations should a priori be solved in a complete single-particle basis while at the HF 
level only occupied states are necessary. In practice, specific methods should be used to truncate 
the basis. A second additional problem is the effective force that should be used in the pairing 
channel. Though zero range forces are clearly very useful at the HF level, they lead to ultra-violet 
divergences in the pairing channel. In practice, either a finite-rang e interaction has to be used or a 
specific regularization or renormalization scheme should be applied Dob84l . lBul02al . rBul02bl . Bul03l | . 



Although these methods provide reasonable solutions for the static case, their application to nuclear 
dynamics is not straightforward. 

Nowadays, the TDHFB theory has essentially b een applied in the small amplitude limit (lead- 
ing to the so-call ed QRPA theory RinSOl . lKha02l | ) or semi-classical limit (hydrodynamical mod- 



interaction 



els TorSTl . lAbrOG 1 ) . Only recently, the dynamical problem has been solved using the nuclear Gogny 



HasO 



B. When is the independent particle approximation valid ? 

In previous sections, we have considered dynamical evolutions where the many-body state is 
forced to remain in a specific class of trial state (Slater determinants or more general quasi-particle 
vacua). This assumption could only give an approximation of the exact dynamics and we do 
expect in general that the system will deviate from this simple state hypothesis. We give here 
some motivations of the introduction of theories beyond mean-field. 

Let us recall the original goal which is to describe as best as possible the dynamics of a self- 
bound complex quantum system. We first assume that a system is initially properly described 
by a Slater determinant"^", i.e. \ip{t = 0)) = with = U^^^aa |— ) where the label a refers 
to initially occupied states (hole states). The exact dynamics of the system is given by the time- 
dependent Schroedinger equation (Eq. ([5])). Unless the Hamiltonian contains one-body operators 
only, the mean-field theory can only approximate the exact evolution of the system. 



1. Decomposition of the Hamiltonian on particle-hole (p-h) basis 

To precise the missing part, we complete the occupied states by a set (possibly infinite) of 
unoccupied single-particle states (also called particle states) labeled by a (associated to the cre- 
ation/annihilation dg and d^). The completed basis verifies 

^|a)(a| + ^|Q)(Q| =/>+(l-/>) = i. (82) 

a a 

From the above closure relation, any creation operator associated to a single-particle state \i) 
decomposes as 

a a 

The particle-hole basis is particularly suited to express any operator applied to the state |(/>) due 
to the properties 

al\<^) = d«|0) = 0. (84) 



The discussion below can easily be generalized to an initial quasi-particle vacuum 
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For instance, restarting from the general expression of H (Eq. ([5])) the different single-particle states 
{i,j,k,l) can be expressed in the particle-hole basis (Eq. (jSSp ). Then, using anti-commutation 
relations dMl) et ^M^, Vijki = -Vjiki = -Vijik = vjuk and Tri2 (ui2/9iy02(l - -P12)) = 2Tri2(-ui2/0iP2) 
as well as (I20D, dll]), ([35]) et dUD, we finally end with 

H\(t)) = \E[p] + Y,-^ h[p]aaalaa HmfIp] 



+ lY.af3af3 Va0f3a aWpf^paa Vres[p] ^^^^ 



where E[p] = {(f)\HMF\4') = Tr [p (t -|- ^?7[/o])] corresponds to the Hartree-Fock energy. 

The above expression is helpful to understand the approximation made at the mean-field level. 
In previous section, we have shown that mean-field provides the best approximation for one-body 
degrees of freedom using the Ehrenfest theorem. In append ix [Gl we show that the mean-field 



evolution can be deduced using the Thouless theorem [Tho6ll | and an effective Hamiltonian where 
Vres is neglected. 



2. Limitation of the mean-field theory 

In expression ()85p . a clear separation is made between what is properly treated at the mean- 
field level {Eq[p] and Hmf[p]) and what is neglected, i.e. V^es[/o]- The latter is called 'residual 
interaction'. At this point several comments are in order: 

• The validity of the mean-field approximation depends on the intensity of the residual inter- 
action which itself depends on the state |^>) and ther efore w ill significantly depend on the 
physical situation. Starting from simple arguments f Lic76l | . the time tsd over which the 
Slater determinant picture breaks down could be expressed as: 

= ^(^ E l(«^l^l«/5>l')~'^'- (86) 

In the nuclear physics context, typical values of the residual interaction leads to tsd — 
100 — 200 fm/c. Therefore, even if the starting point is given by an independent particle 
wave-packet, the exact evolution will deviate rather quickly from the mean- field dynamics. 
This gives strong arguments in favor of theories beyond TDHF in nuclear physics. 

• An alternative expression of the residual interaction which is valid in any basis, is 

Vres[p]l2 = ^(1 - - P2)l^l2PlP2- (87) 

This expression illustrates that the residual interaction associated to a Slater determinants 
could be seen as a "dressed" interaction which properly accounts for the Pauli principle. 
Physically, the residual interaction corresponds to direct nucleon-nucleon collisions between 
occupied states (2 holes) which could only scatter toward unoccupied states (2 particles) due 
to Pauli blocking. We say sometimes that the residual interaction has a 2 particles-2 holes 
(2p-2h) nature. 

Due to the residual interaction, the exact many-body state will decompose in a more and more 
complex superposition of Slater determinants during the time evolution. As stressed in the in- 
troduction of this lecture, due to the complexity of the nuclear many-body problem, the exact 
dynamics is rarely accessible. In the following section, methods to include correlations beyond 
mean-field, like direct nucleon-nucleon collisions or pairing, are discussed. 
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C. General correlated dynamics: the BBGKY hierarchy 



Using the Ehrenfest theorem (section III A4p , we have shown that the mean-field theory is par- 
ticularly suited to describe one-body degrees of freedom. A natural extension of mean-field consists 
in following explicitly two-body degrees of freedom. Considering now the Ehrenfest theorem for 
the one and two-body degrees of freedom leads to two coupled equations for the one and two-body 



density matrix components pi^ = {a\ai) and p. 



it; 



,(2) 



^^^Pl2 = [tl +t2 + \vi2,pi2\ + ^Trs [Wl3 + V23, PUs] ■ 



(88) 



Above equations are the two first equation of a hi erarchy of equations, known as the Bogolyubov- 
Born-Green-Kirkwood-Yvon (BBGKY) hierarchy Bog46l . Bor46l . Kir46l | where the three-bocly den- 
sity evolution is also coupled to the four body density and so on and so forth. Here, we will restrict 
to the equations on p^^^ and p^^-' which have ofteri served as the starting point to develop transport 
theories beyond mean- field [CasQd . lRei94l . lAbe96l . lLac04i | . 



D. The Time-Dependent Density-Matrix Theory 

Previously, we have shown that the mean-field dynamics neglect two-body and higher correla- 
tions. Then the equations on p^^^ reduce to TDHF. A natural extension which includes two-body 
effects is to treat explicitly two-body correlations and neglect only three-body (C123 = 0) and 
higher correlations^^. The resulting theory where the one-body density pi and the two-body corre- 
lation C12 are followed in time is generally called Time-Dependent Density-Matrix (TDDM) theory 
(see for instance [Cas90]) 



ih^^pi 



[hi[p], pi] + ^TT2 [vi2,Ci2] 



12 



[hi[p] + h2[p],Ci2 



Pl){l - P2)vi2PlP2 - PlP2Vl2{l - Pl)(l - 
Pi - P2) ^12Cl2 - Ci2-Ul2 (1 - Pl - P2) } 



92] 



+Tr3 [via, (1 
+IV3 [7)23,(1 



^13)^1^23(1 

i^23)piC23(l 



P12)] 
P12)] . 



B12 
Pl2 

H12 



(91) 



where we have dissociated explicitly three terms which will be responsible for the build up of 
correlations in time. The Born term, B12, contains the physics of direct in-medium nucleon- 
nucleon collisions. Comparing B12 and Eq. ()87p . we see that it is directly proportional to the 
residual interaction. Indeed, starting from a Slater determinant (Ci2(to) = 0), this is the only 
term that do not cancel out in the evolution of C12 over a short time scale. In particular, it will 
be responsible for the departure from an independent particle picture. The physical interpretation 



Introducing the permutation operator P12 between two particles, defined as P12 = The two-body correla- 
tion matrix is given by: 



C12 = P12 - pip-zil - P12) 
while the three-body correlations C123 reads 

C123 ~ pi23 — P1C23 (1 — P12 — Pis) — P2C1S (1 — P21 — P23) 

-P3C12 (1 - P31 - P32) - P1P2PS (1 - P13) (1 - P12 - P23) ■ 



(89) 
(90) 
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of the term P12 and H12 is less str aightfor ward. For instance, it has been shown that P12 could be 
connected to pairing correlations Toh04| (see discussion below) while H12 contains higher order 
p — p and h — h correlations. It is finally worth mentioning, that th e last t erm could eventually be 
modified to better account for conservation laws (see discussion in Pet94l ]). 

Application of the TDDM theory faces two major difficulties. First, since we are consider- 
ing explicitly two-body degrees of freedom, we have to deal numerically with huge matrices and 
appropriate truncation schemes should be performed. Second, numerical applications are mainly 
possible with contact interaction (Skyrme like). These interactions, which are zero range in r- 
space are thus of infinite range in momentum space. This unphysical behavior of the interaction 
is critical in practice, since during nucleon-nucleon collisions, particles may scatter tow ards to o 
high momentum. No clear solution to this problem exists so far in the TDDM theory Lac04 1 . 
Due to these difficulties, only fe w applications have been carried out so f ar for col lective vibrations 
Deb92l . lLuo99l . iTohOll . lToh02a | , and very recently for nuclear collisions Toh02bl | . 



E. Link between TDDM and TDHFB 



The connection between the TDDM and TDHFB has been clarified in Ref. TohOJ]. Assuming 
a separable correlation in the p-p and h-h channels given by equation (|71|) gives 



^ (A |Tr2 [vi2, Ci2]\ A') = ^ (AA; \vi2\ mn) (mn \Ci2\ X'k) - {\k \Ci2\ mn) (mn \vi2\ \'k) 



kmn 



XX' 



(92) 



where A is nothing but the pairing field introduced above. Then, the one-body density evolution 
reduces to 



ih^p = [h[p],p] + kA* - Ak* 
at 



(93) 



In Ref. TohOJ], it has been shown that neglecting B and H in the second equation of (|9ip leads 
to a TDHF like equation. Keeping only P and assuming (j7ip leads to 



d 



.4- ^ ■ J- ( djKij ^ d^hl 1 

dt ^ I dt ^ dt / 



mnpq 

^ ' ^ j i^ml^mk "I" ^km^ml) + ^ ^ l^mn'^mnpqi^kp^lq ~ ^kpPql ~ ^IqPpk 



mnpq 



(94) 



where we have identified the terms proportional to Kij et k^^. We then get 

ih OfKij — ^ ^^ {hjm^mj ~l~ ^im^mj) ~^ 7^ ^ ^^ ij^ijmn PipVpjmn Pjp^ipmn^ l^mn^ 



(95) 



Using the expression of the pairing field (Eq. (j76p ). we finally recover the TDHFB equation (j79p . 
The above equation does not insure that the correlation matrix remains separable during the time- 
evolution. However, assuming that Cijki{t) ~ Kij{t)K*^i{t) is valid for all time, the equations of 
motion identify with the TDHFB equation. It is worth mentioning that the above technique gives 
an alternative derivation of the TDHFB equation starting from TDDM which in addition illustrates 
the physical content of P. 
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F. Extended and Stochastic Time-Dependent Hartree-Fock 

The pairing correlations become less important when the internal excitation of the system 
increases. In this case, direct nucleon-nucleon collisions are expected to dominate correlations 
beyond mean-field. A possible way to treat the latter and to escape part of the complexity of 
TDDM is to focus on the one-body density while the effect of correlations on this quantity is only 
treated approx i matively. This could be done by neglecting the terms P and H in the evolution 
of Ci2 [Won78l . IWonTd . lDan84 lBot9Cll . I Avi8Cll | . Then, the equation for the two-body correlations 



takes a simple form 

i^Q-^Cu - [hi[p] + h2[p], Cu] = Bi2. (96) 

Solving this equation formally, we can develop the correlations over a time interval from an initial 
time to to a time t as 

Ci2{t) = -'- [ds Ui2 {t, s) Si2 (s) Ul^ (t, s) + 6Ci2{t) (97) 
Jto 

where U12 = Ui®U2 represents the independent particle propagation of two particles with U{t, s) = 
exp ^"■^ j h[p{^')\dt'^ ■ In expression ()97p . the first term in the right hand side represents corre- 
lations developed by the residual interaction during the time interval. The second term describes 

propagation of the initial correlations Ci2(to) from to to t, i.e. dCi2{t) = ^712(^5 io)C'i2(io)f^i2(^' *o)- 
Reporting this expression in the evolution of pi, we end up with a generalization of the TDHF 
theory (where we omit the label "1" in pi) 

ih^^p=[h[p\,p]+K[p\ + 6K{t) (98) 
where K[p\, called collision term, reads 

K[p] = [ds TT2[vu,Ui2{t,s)Bi2{s)ul^{t,s)] (99) 

while 5K{t) is given by 

5K{t) = Tr2bi2,5Ci2(t)]. (100) 

The term 6K{t), which accounts for the initial correlation Ci2(to), dominates around to. This 
term, in principle, contains all order correlations that are accumulated up to time to, and it is a 
priori a very complicated quantity. Its treatment is clearly out of the scope of a one-body transport 
theory and statistical assumption is generally made on the initial correlations. It is assumed that 
the exact two-body correlations accumulated until to exhibits random fluctuations. As a result, 
the average value of the initial correlations vanishes. This assumption is known as the "molecular 
chaos assumption" in classical transport theory and it c o rrespon ds to factorization of two-particle 
phase-space density before each binary collision Kad62l . lHua62l l . Eq. ([US]) is then replaced by an 
ensemble of one-body evolutions 

ih-^-^p- = [hnp^]+K[p^]+5K"it) (101) 

where "n" now refers to the particular stochastic path and where 6K'^{t) is a fluctuating oper- 
ator which vanishes in average. The complete statistical description of 5K'^{t) can be found in 
Refs. AviOll . lLac04i |. Eq. (jlOip is the starting point of most of the microscopic transport theories 



that are applied nowadays in Heavy-Ion reactions at intermediate energies. However, due to its 
complexity, mainly the semi-classical version of Eq. (jlOip . known as Bolt zmann Langevin theory, 
has been applied to realistic situations (for a recent review, see |Cho04| l. Such a semi-classical 



approximation is not expected to apply in low energy nuclear reactions and we will concentrate 
here on its quantum version. 
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1. Average evolution: irreversible process in Extended TDHF 



We first concentrate on the average evolution and illustrate advantages of the introduction of 
a collision term on top of the mean-field dynamics. Averaging Eq. (I98p over different trajectories 
leads to an equation where is replaced by the average one-body density, denoted by p and 
where only the collision term K[p] remains. The resulting theory is called Extended TDHF with a 
non-Markovian collision term (or with "memory effects"). The terminology " non-Markovian" (in 
opposition to " Markovian" ) comes from the fact that the system at time t depends not only on 
the density at time t but also on its full history due to the presence of a time integral in Eq. ([UU]) . 

At any time, eigenstates oi'p(t), denoted by \a{t)) could be found. In this basis, called hereafter 
natural basis or canonical basis, p(t) reads 



p{t) = Y^Ht)) n^t) {a{t)\ 



(102) 



Using the weak coupling approximation in combination with the first order perturbation theory, the 
ETDHF equation can be transformed into a generalized master equation for occupation numbers 
which accounts for the Pauli principle 



di 



n, 



t 



to 



{ (1 - n„ (s)) W+ (t, s) - n„ (s) W" (t, s)} . 



(103) 



where the explicit form of the gain and loss kernels can be found in Ref. Lac99l | . There- 
fore, in contrast to TDHF where occupation numbers are constant during the time evolution, in 
ETDHF the evolve and could eventually relax toward equilibrium. Such a relaxation is the only 
way to properly account for thermalization process in nuclei. The inclusion of correlation effect 
with Extended TDHF has been tested in the simplified case of two interacting nucleons in one di- 
mension |Lac99l |. In this case, the exact dynamics could be solved numerically. In Fig. [T5l starting 
from an initially uncorrelated state, the exact evolution of single-particle occupation numbers is 
compared to the Extended TDHF prediction. Fig. [T5l shows that the Extended TDHF is able 
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Figure 15: Exact (circles), TDHF (dashed line) and Extended TDHF (solid line) evoluti ons of a two nucleon 
correlated system. Left: Occupation numbers, right: one-body centroid (adapted from [Lac99j . 

to reproduce fairly well the evolution of one-body occupation numbers and one-body observables 
over long time. This result is very promising and indicates that Extended TDHF seems to be an 
appropriate theory for the description of dissipation when the residual interaction is weak. The 
above application has also demonstrated the importance of memory effects to properly describe 
quantum interacting systems although it significantly increases numerical efforts. 
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2. Discussion on Stochastic mean-field dynamics 



Due to the underlying hypothesis of Extended TDHF, such an average theory is expected to 
apply only over short time scale. The long time dynamics requires to include both dissipation and 
fluctuations and therefore to solve explicitly the stochastic version of the transport theory (Eq. 
(jlOip ). Up to now, stochastic mean- field has been essentially applied in the semi-cl a.ssical limit 



AviOll ] 



was 



neglecting explicitly non-Mar kovian effects. Only recently, the theory proposed in Re f. ^ ^ 

used to study small amplitude collective vibrations in quantum system [L acOl , lLac04l | including full 
memory effects. In that case, the description of damping of giant resonances is greatly improved by 
including both dissipative and fluctuating kernels. However, due to numerical as well as conceptual 
difficulties, the application of the above stochastic quantum mechanics to large amplitude collective 
motion, like nuclear collisions, remains an open problem. 

At this point, it is worth mentioning t hat exte nsive works have been devoted to the formal deriva- 
tion of dissipative quantum mechanics Kad62l | a nd/or r e lated s t ochast i c equat i ons for f ermions , 
including Markovian and non- Markovian effects [AviSOl . iBalSll . lAviSSl . lAviOll . lRei92al . lRei92bl . 
Ohn95| . These approaches have in common that the residual part of the interaction introduces 
disorder on top of the mean field. Such theories end with rather complex transport equations 
which are hardly applicable in realistic situations like Heavy-Ion collisions. Connections have been 
made recently between the dissipative equations of the one-body density and quantum jumps be- 
tween Slater de terminants using either p erturba tion theory or the connection between ETDHF 
Rei92al . lRei92bll and Lindblad equations Lac06al ] generally found in the theory of open quantum 



systems Bre02l | . These theories, in which the many-body density are explicitly replaced by an av- 
erage over Slater determinant densities, i.e. D = {(j)"') (</!)"■ | where each Slater determinant evolves 
according to a Stochastic TDHF equation, open new perspectives. 



G. Functional integrals and exact treatment of Many-Body correlated systems with 

Stochastic mean-field theories 

1. General Discussion 



Starting from a simple Slater determinant state \ip{t = 0)) = \(f){to)), correlations will develop 
in time and we do expect that the exact Many-Body state writes 



(104) 



where \(f)k) denotes a complete (eventually time-dependent) basis of Slater-Determinant states. 
Accordingly, the many-body density writes (see appendix IC]) 



m = Y,ckit)ci,it)\Mt)){^k'it)\. 



(105) 



k,k' 



The extended and stochastic version of TDHF presented in sections IIVF[ implic itly assume that 
the many-body density can be properly approximated by its diagonal components Rei92al . lLac06a( | 



Dit) ^Pfc|(/>fc(t))(</.fc(i)| 



(106) 



where Pk = |cfc(t)p. Then a probability Pk obeys master equations which eventually could be 
simulated by quantum jumps. The resulting density is obtained through the average over different 
stochastic paths, i.e. 



D{t) ^ \Mt)){Mt)\- 



(107) 
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Physically, this could be understood as follows. The irrelevant degrees of freedom (complex internal 
degrees of freedom) interact with the relevant degrees of freedom (single-particle degrees of freedom) 
and induce a fast decay towards ze ro of the off-d iagonal matrix elements. This phenomenon 



is known as a decoherence process jKue73l . lKie03l ] . It is clear, that such an approximation is 



associated with a loss of some quantum effects such as interferences between different channels. 
Therefore we do expect that most of the extensions of TDHF presented above, as TDHF, will 
miserably fail to account for most of the true Many-Body quantum effects. 

The goal of the present section is to demonstrate that one could always treat exactly the 
Many-Body density given by Eq. (jl05p using an appropriate stochastic process between Slater 
determinants. Since there is no free lunch, the exact density is then obtained by an average 



D{t) ^ \Mt))mt)\ (108) 
where states in the left differ from states on the right. 



2. Functional Integral technique in a model case 
Functional integrals techniques have often been used to repla.ce the exact Ma ny-Body problem 



l^unctional mtegrals tecUmques nave oiten been used to repia.ce tne exact Ma ny-.l3ody problem 
by an average over different " effective" one-body problem LevSOal . [LevSObl . [Neg88l | . In Ref. jKoo97 |. 



the general strategy to obtain ground state properties of a many-body system using Monte-Carlo 
methods, the so called Shell-Model Monte-Carlo, is described. Recently, this technique has been 
combined with mea n-field t heory to obtain Stochastic TDHF equations which in average leads to 



the exact evolution CarOll . |Jui02l | . We give here a brief description of the method. 



3. Simple Introduction to functional integrals 

We again consider that at a given time, the Many-Body state is a Slater Determinant \^{t)) = 
{(p). For a short time step At, we have 

m + At)) = exp (^^H^ \cPit)) ^(l + ^H + o(At)) \<p{t)) . (109) 

Due to the presence of a two-body interaction in H, the state \ip{t + At)) differs from a Slater 
Determinant. Let us show how it can be written as a sum of Slater Determinants. For simplicity, 
we first assume that H = Hi + O^, where Hi and O are both one-body operators. Therefore 
corresponds to a two-body interaction. We introduce the notation G{x) for a normalized gaussian 
probability of the variable x with mean zero and variance 1 

dx xG{x) = 0, x2 = / dx x^G{x) = 1. (110) 

-oo J —oo 

We define the complex number Aw = \J'^^ and the one-body operator ^(At, x) written as 

S{At,x) = ^Hi + xAujd. (Ill) 
in 

Considering the average value of 5(At, x) and keeping only terms up to At, we obtain 

/ + 00 , A J. 1 A J. 

dx e^(^*'^)G(x) = 1 + =^Hi + X AcoO + -x^ {Au;)^d^ + o(At) = 1 + =^H + o(At). (112) 
in 2 in 
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By averaging over the different realizations of x, we recover the exact propagator over a short time 
step. Note that more general relations could be found using the Hubbard-Stratonovish transfor- 
mation (see for instance [Koo97]). Using the above relation, we see that 

—H\\^(t))= dx G(x)e^(^*'^)|0(t)) = / dx G(x)|</.^.(t + At)). (113) 

/ o — OO "/ —(X) 

Due to the one-body nature o f S, each \4)x{t + At)) is a Slater determinant. Indeed, according to 
the Thouless theorem Tho6ll | (see appendix [Gl an exponential of a one body operator transforms 



a Slater determinant into another Slater determinant. We have therefore demonstrated that the 
evolution of the exact state could be replaced by an ensemble of Slater determinants. The method 
could then be iterated for each \(j)x{t + At)) to obtain the long-time evolution as a superposition 
of independent particle states. 
Several comments are in order: 

• Since S'(At, x) is not a priori Hermitian, the dynamics does not preserves the orthogonality 
of the single-particle wave -funct i on. Su ch a non-orthogonality should properly be treated 
during the time evolution Jui02l . iLacOSl ] . 



Starting from a Many-Body density written as D{t) = \(j)) at an intermediate time, the 
average density writes 



D{t) = \Ut)){Ut)\ (114) 
with two states, where |0i) evolves according to Eq. (jllSp while (021 evolves according to 

(</.2(t + At)| = (02(t)|exp|-^iJi+yAu;*o| (115) 

where y is a noise independent of x, with mean zero and yy = 1. Since the evolution is exact, 
any one-, two- or k-body observable A estimated through (A) = Ti{D(t)A) will follow the 
exact dynamics. 



4- General two-body Hamiltonian 

In previous section, we have illustrated how the exact dynamics of a system could be replaced 
by a set of stochastic evolutions of several Slater determinants using a schematic two-body Hamil- 
tonian. We now start from the general expression of the Hamiltonian (|85p to introduce the exact 
Stochastic Mean-field (SMF) technique. In this expression, the Hamiltonian is naturally splitted 
into a mean- field part Hmf and a two-body residual interaction V^es[/o]- It tu rns out that this 



interaction can always be decomposed as a sum of square of one-body operators jKoo97l | 



(a^|^i2|a/3) = J^CA(a|dA|a)(^|OA|/?), (116) 
A 

where Oa is a one-body operator while c\ denotes a set of constants (eventually complex). The 
discussion in the previous section can then easily be generalized. The residual interaction Vres 
factorizes as 



where Oa = Ylaa('^\^^\'^)^e'^a- Therefore, for realistic interactions one should introduce as many 
stochastic Gaussian independent variables as the number of operators entering in the sum. In 
practice, this number defines the numerical effort which in general is very large. For this reason, 
only few applications to the dynamics of rather simple systems exist so far. Finally, the extension 



of above stochastic theories to HFB states can be found in Ref. [Lac06bl | while an explicit link with 
observables evolutions is studied in Refs. LacOTl . [LacOSl ] . 
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H. Summary 

In this section, wc have summarized some of the possible ways to extend TDHF, some of them 
are able to incorporate pairing correlations (like TDHFB or TDDM) while others concentrate 
on direct nucleon-nucleon collisions (ETDHF). Though promising, these theories to the nuclear 
many-body problem remain very challenging to apply. 

A second difficulty which has been largely hidden in the present lecture is that all applications of 
dynamical quantum transport theories, like TDflF, to nuclear reactions was possible only because 
of the introduction of effective interactions (essentially Skyrme like). These interactions have led 
to the more general concept of Energy Density Functional (EDF) and are expected, in a similar 
way as the Density Functional Theory (DFT) in condensed matter, to incorporate most of the 
correlations already at the mean- field level. Then, the notion of "beyond mean- field" calculations 
becomes ill defined. For instance, all theories (extended, stochastic, functional...) start from a 
Many-Body Hamiltonian. In the EDF context, such an Hamiltonian, although it exists, is not 
simply connected to the EDF itself. As a consequence, the Hamiltonian derivation could serve 
as a guideline but a proper formulation in the EDF framework is mandatory. Large debates 
exist nowadays on the validity and foundation of the nuclear EDF applied to static properties of 
nuclei. Similar discussion should be made on what we should definitively call Time-Dependent 
EDF (TDEDF) and not "TDHF". 
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Appendix A: BASICS OF QUANTUM MECHANICS 

Here, we summarize some aspects of quantum mechanics and more specifically second quanti- 
zation. Only concepts useful for the the present lecture notes are introduced (see jMes59l . lBla86l | 
for more details). 



1. Single-particle creation/annihilation operators 

The single-particle wave- function with spin and isospin components denoted respively by s 
and r can be written as 

^r(r) = {rsT\i), (Al) 
The associated creation operator a\ acting on the vacuum |— ) verifies 

al|-) = |i). (A2) 

Its Hermitian conjugated is the annihilation operator ai\i) = |— ). For all annihilation operator, we 
have 

a,|-) = Vi. (A3) 

This property defines the particle vacuum associated to the set of operators and a. 

We consider here fermions. The Pauli principle imposes that two identical fermions could not 
be created in the same state a}^\i) = 0. Creation/annihilation operators verifies anti-commutation 
rules 

{a|,d]} = {a„aj} = (A4) 
{a|,a,} = % (A5) 

where {a, h} = ab + ha. 

It could be useful to go from one single-particle basis to another, and to express cre- 
ation/annihilation of the first basis in terms of creation/annihilation of the second basis {|^)}. 
To do so, the two basis should obviously represent the same Hilbert space (and be complete with 
respect to this space). Fort instance, the relation of any single-particle basis to the coordinate 
representation is 

airsT) 



5^^r(r)a, (A6) 

i 

a\vsT) = Y,^r{r)a\. (A7) 

i 

Using the closure relation X^st/'^^ |rsr)(rsr| = 1, previous relations can be inverted as 

ai = ^ /dr (^r*(r) a{vsT) (A8) 
al = /dr V^r(r) at(rsr). (A9) 
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2. N identical particles state 

We denote by l^/') a A^-body state which is a priori correlated. Occupation numbers of a 
correlated state, in opposite to an uncorrelated state, have values which might differ from zero 
and one. A A^ independent particles state corresponds to a specific case where A" occupation 
numbers are exactly one while the other are zero. We not such a state It can be written as an 
antisymmetric product of single-particle states (called Slater determinant due to its specific form, 
see Eq. (RMll ) 

= \<P,,...,^) = i |1 : z/i, 2 : N:i^n), (AlO) 

where |1 : a, 2 : /3, • • •) means that the particle 1 is in the state \a), particle 2 occupies state 
The "v/aT coefficient insures the proper normalization of the wave while the anti-symmetrization 
operator writes 

A = ^ E sign(P)P. (All) 
permutation P 

Here sign(P) = 1 (resp. -1) for odd (resp. even) permutations of particles. For instance, its action 
on a two-particle wave-function reads 

^ |1 : Q,2 : /?) = ^ (|1 : Q,2 : /3) - |1 : /3,2 : a)) . (A12) 

In second quantization, an independent particle state reads 




(A13) 



3. Wick's Theorem 

The Wick's theorem enables to simply express expectation values of observables for specific 
states (that could be considered as vacua). Let us first define the term contraction. Consider two 
creation and/or annihilation operators, denoted by A and B (these operators could also be linear 
combinations of creation/annihilation operators) and |0) the associated vacuum. In practice, this 
vacuum could be a single-particle vacuum, a HF state denoted by \(p) or a quasi-particle (HFB) 

vacuum (see section IIV Ap . The contraction AB is defined as the expectation value of AB on the 
vacuum 

Ab = {0\Ab\0). (A14) 

Then, the Wick's theorem states: 

The expectation value of a product of creation/annihilation operators on their vacuum |0) is 
equal to the sum of all possible products of contractions of pairs of operators, each product of 
contractions being multiplied by + or - depending on the parity of the permutation needed to put 
together the contracted operators. 

Let us illustrate this on the overlap between two anti-symmetric two independent particles states 

{<Pai3\<Pfj.u) = {-\df3aaal,al\-). (A15) 



Contractions are given by ajOj = ajcij = and ajoj = {i\j). We then deduce 



{(paisl^Pfiiy) = aisdl Qaol - CLpal aaol = {/3\u) (a|/Li) - (/3|/u) {a\u). (A16) 
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This result can be generalized to the overlap between two A^-particles Slater determinants 
constructed from two different basis and {|a)} 



l---ajv) — ( 



(A17) 



We can then use the Wick's theorem to express this overlap in terms of a determinant made of 
overlaps between single-particle states (z^jlaj) 



(A18) 



4. Basis of N particles states 

Starting from a single-particle complete basis, we can always construct a basis of A'^-body Slater 
determinants. For instance, we calculate the overlap between two Slaters formed from a priori 
different single-particle states of the same basis. According to Eq. (jAlSp . we have 



(A19) 



This overlap is ±1 if 



and 



contain exactly the same occupied states and if at 



least one of the single-particle states is different. Therefore, two different Slater determinants are 
orthogonal. It could be shown finally that the set of Slaters constructed by picking up all possible 
combination of N single-particle states forms a complete basis of the A^-body Hilbert space. In 
other words, any many-body state, correlated or not, can be written in the form 



(A20) 



5. Independent particles wave- functions 

The Wick's theorem is also helpful to express the wave function of a A^ particles state. We want 
to write a wave function in position-spin-isospin representation. Associated single-particle states 
are = |r„s„T„), where r is the position, s and r the spin and isospin projections respectively. 
Using Eq. (lAln]) . a Slater built using this basis reads 



\i^...i^) = ^\A\\:ix---N:v^). 
The A^ particles wave function is then given by 



(A21) 



N 



For an independent particle system, the wave function reduces to 

1 



(-|a(^Ar)---a(6) al ■■■al \-) 



(A22) 



(A23) 
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From Eq. (jAlSp . we see that we simply obtain the standard Slater determinant formula 



(A24) 



Let us check that this wave- function is properly normalized. Introducing the notation Jd^ = 
YlisT I^^-i nomalization reads 



Norm[</>] = jdii---diN (Ptj_...UNi^^ ' " ^n) (pu^.M^iCi ' ' ' Cn) 



Nl 



d6---de 



N 



(A25) 
(A26) 



Developing the determinants and using Jd^ 'fiiiO VjiO = ^ji '^^ finally get 

Norm[(/>] = ^ JdCi---dCN Yl PWIACi)^uACi)---V>:A^n)v>uACn)} 

Permutation P 

= 1. (A27) 



Finally, it is worth to mention, that, using Eq. (|A22p and (]A25p . the closure relation for the 
particles states reads 



iiv = ^ ldCi---d^N \^i---Cn){Ci---^n\. 



(A28) 



Appendix B: ONE-BODY OBSERVABLES 

One-body observables are generally written as a sum of operators acting on each single-particle 
components i independently, i.e. 



N 



F = Yhi). 



i=l 



In second quantization, one-body operators reads 



j) aldj. 



(Bl) 



(B2) 



Both expressions lead to the same action on a A^-body wave function. For instance, considering 
the independent particle case, using Eqs. ()Bip . (jAlOp . ()A12p and the notation f\j) = fij\i), we 
have 



TV 



i=l 
N 



i=l 
N 



Y Yl V, «li • • • "Li «M «Ui • • • "^iv 



(B3) 



4=1 /X 
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Eq. dHl]) leads to 

F\(t>u^-ut,) = ^ ff,ual al^--- alj-)- (B4) 
Only occupied u contribute due to Eq. IA3I and we get 

N 
N 

= ft^'^^ ■ ■ ■ 4 «Ui • • • «L I-) (B5) 

where we have also used Eqs. ()A4p and (jASp . This result is nothing but Eq. (|B3p . The table [U 
gives few examples of one-body operators commonly used in the nuclear reaction context. 



Observable 


Standard form 


Second quantization 


Center of mass position R 




^Esr /dr r a^(rsr) a(rsr) 


Center of mass momcntumP 




Esr/dP P a^psT) a{psT) 


Particle number N 






Monopole operator Qo 




^Esr/di- a^rsT) a{rsT) 


Quadrupole operator Q20 




\f^J2srI'^^ {2z-x-y) at(rsr)a(rsr) 



Table II: Typical example of one-body operators. 



Appendix C: DENSITY MATRICES 

The A^-body density matrix of a A^ particle system (described by {ip)) contains all the infor- 
mation on the system. It could be written as an operator acting on the A^-body Hilbert space 
as D = \ilj){ip\ (see for example |Sur95l . lAbe96l . lLac04| ). If we are only interested in the M-body 



observables (M < A^), we do not need all the information carried by D. A natural way to reduce 
the information is to introduce the M-body density matrix of the state It could be defined 
from D as 

Equivalently, the M-body density matrix components are defined by the relation 
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In coordinate space = {rsr}), it writes 



{N -M)l 

Other notations, like p^'^'^^ = pi,,,M = p{^---M), are also sometimes used. The main advantage of 
the M-body density is that it contains all the information on the system as far as M-body or lower 
observables are concerned. For instance, the expectation value of a M-body observable O^^^ reads 

(oW)^ = Jdc,...dCMde^...deM p^^'\ei-eM,^i-M oW(6...eM,ei...c;^) 

= TVi,..M[p(*-')0(^)]. (C4) 

Appendix D: TWO-BODY CORRELATIONS 

When we want to get the expectation value of a two-body observable, as it is the case for 
the Hamiltonian H, we do a priori need the two-body density matrix p^^^ of the system (see 
appendix E]) . The two-body density matrix can always be decomposed into an uncorrelated part 
(anti-symmetric product of one-body densities) and a correlation part, denoted by C^'^^ 

Pijki - ^k'^^'^j)^ - Pji Pik Pa Pjk+'-'ijki \^^) 

The correlation operator C*^^-* can be seen as the part of the two-body density which cannot be 
written as a product of single-particle operators. It is also possible to use the notation 

P12 = P1P2 (1 - P12) + C12 (D2) 

where P12 corresponds to the permutation operator between particles 1 and 2. We recall that 
notations O*-^-* = 0(1,2) = O12 are equivalent. Note that higher M-body correlations matrices 



(with M > 3) could be defined in a similar way [LacOj]. 



Appendix E: HARTREE-FOCK AND QUASI-PARTICLE VACUA 

A Slater determinant is sometimes called HF vacuum. The state |0) is indeed a vacuum for 
the creation operators bji and annihilation operators 6^ written, in the basis that has served to 
construct the many-body state, as 

^Jt = (1 - n-f,) al + (El) 
6^ = (1 - n^) + rif, dj^ (E2) 

where = 1 for occupied states (also called hole states) and for unoccupied states (particle 
states). It could easily be checked that these new operators verify Eqs. ()A4p and (|A5p and that 
10) is a vacuum for them, i.e. 

b^\<P) = V/i. (E3) 

Slater determinants are a specific case of quasi-particle vacua. More generally, quasi-particle 
creation/annihilation operators are defined through a linear combination of single-particle cre- 
ation/annihilation operators (a|,aj) 

= EiU*aat + V*^dl ^^^^ 
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where matrices U etV are chosen in such a way that the quasiparticle operators verify the fermionic 
anti-commutation rules (Eqs. ()A4IIA5P ). The vacua associated to these operators can be written 
as 

li^bogo) =Cl[Pk\-) (E5) 

k 

where C is a normahzation constant. We see that the latter expression insures $i\ipbogo) = 0. 



Appendix F: TWO-BODY DENSITY FOR INDEPENDENT PARTICLE SYSTEMS 



We show here, using the Wick's theorem, that it is equivalent to have the property C^'^^ = 
and to consider a system of independent particles. Starting from a Slater determinant |(/>), using 
the fact that this state is a vacuum for the 6^ (see appendix [E|), and the inverse of Eqs. (|EmE2p 
in a given basis denoted by \i) 



4 



^ {n\i) bf, + {1 - Uf,) bl 
^ Uf, bl + {1 - Uf,) b^ 



(Fl) 
(F2) 



with = or 1, we can use Wick's theorem to express two-body matrix elements 



(F3) 



Here contractions are made with the HF state \4>). Using a|aj = {a\aj)^ = pji and a|aj = 0, we 
finally obtain that the two-body matrix elements read 



.(2) _ .(1) .(1) _ .(1) .(1) 
rijkl ~ Pjl Pik Pil Pjk ■ 



(F4) 



Here, the two-body correlation C^"^^ is strictly zero. In the notes, we often used contractions in the 
specific basis of occupied states, which simply read 



(F5) 



Appendix G: MEAN-FIELD DYNAMICS FROM THE THOULESS THEOREM 

Starting from Eq. ()85p and an initial Slater determinant, we neglect the residual interaction, 
i.e. we only keep Hmf[p\ in the dynamical evolution. Then, over an infinitesimal time step dt, the 
new state is approximated by 

\^{to + dt)) ~ exp (^Hmf[p]^ 10) . (Gl) 

where Hmf[p] = E[p] + X]/3o Hplpa ^^^^a, and E[p] = Tr [p (t + |?7[/3])] . We recall that a denotes 
an occupied state and a an unoccupied one. As we will discuss below, to recover the single-particle 
evolutions, it is convenient to rewrite the Mean-Field Hamiltonian in a slightly different manner 
(using d^pCLa \(t>) = Sap 10) ) 

Hmf[p]\<P) = [-^TiipU[p]) + Y,{mp]h)alaa + J2 4«") '"^^ 

= ( - \Tr{pU[p]) + ^ {i\h[p]\a) alda) |0) (G2) 

ia 
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where "i" denotes a complete single-particle basis. The term — ^Tr(pC/[p]) does not influence 
the single particle evolution (we omit its contribution in the following). The propagator ex- 
pressed as an exponential of one-body operators and therefore, according to the Thouless Theorem 



Tho61[, transforms a Slater Determinant into another Slater determinant. Indeed, the fact that 



^-^Hmf^jeHmf = 1 and e^^^/F |_\ |_\ ^llow us to write 



dt Tj ± dt Tj dt Tj J. dt Tj \ dt fj J. dt jj , dt -fj , , 

= eih^MFfj]^^^-JE^MF^Jfi^MFg)^^^-iE^MF . . . g+ Ift -« AfF ^ g" 77-^ A/F g+^i^MF |_^ _ 

Considering the transformation of each creation operator separately, we have 

e^HMFal^-fkHMF = ai + ^[Hmf, ai] + o{dt) 

= + § I^WI + ^(^*) = + ('(dt), (G3) 

i 

where the expression of the mean-field operator defined in Eq. (|85p has been used. From the above 
identity, we see that the propagated many-body state writes \"^{t + dt)) oc nQ,oJj,_(_^^ |— ) where, 
using \i) {i\ = 1, the single-particle states evolve according to 

ih^ = h[p]\a) (G4) 

which is nothing but the standard mean-field evolution. Therefore, we have shown in this appendix 
that the mean-field evolution is recovered by neglecting the two-body residual interaction. 

The fact that we exactly recover the standard single-particle wave-function evolution is due to 
the specific way we have written the Mean-Field Hamiltonian in Eq. ()G2p . If, instead, we take the 
original expression and use the fact that E[p] induces only a phase, then, using the same technique, 
we would get 

ih^ = il-p)h[p] \a). (G5) 



These single-particle equation differs from the standard Eq. ()G4p but contains the same informa- 
tion. Indeed, we have 

ihdtp = ihY, [idt\a)) {a\ + \a) {dt{a\)] = {I - p) h[p\ p-p h[p] (1 - p) = \h[p\rp\ (G6) 

a 



where we have recovered the correct one-body density evolution. 
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